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Abstract: 

We consider some specific inverse problem or "bottom-up" reconstruction strategies at the LHC for both general 
and constrained MSSM parameters, starting from a plausibly limited set of sparticle identification and mass mea- 
surements, using mainly gluino/squark cascade decays, plus eventually the lightest Higgs boson mass. For the three 
naturally separated sectors of: gaugino/Higgsino, squark/slepton, and Higgs parameters, we examine different step- 
by-step algorithms based on rather simple, entirely analytical, inverted relations between masses and basic MSSM 
parameters. This includes also reasonably good approximations of some of the relevant radiative correction calcula- 
tions. We distinguish the constraints obtained for a general MSSM from those obtained with universality assumptions 
in the three different sectors. Our results are compared at different stages with the determination from more standard 
"top-down" fit of models to data, and finally combined into a global determination of all the relevant parameters. Our 
approach gives complementary information to more conventional analysis, and is not restricted to the specific LHC 
measurement specificities. In addition, the bottom-up renormalization group evolution of general MSSM parame- 
ters, being an important ingredient in this framework, is illustrated as a new publicly available option of the MSSM 
spectrum calculation code "SuSpect". 



I. INTRODUCTION 
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If supersymmetry shows up at the LHC, it may be that only a limited part of the predicted Minimal Supersymmetric 
Standard Model (MSSM)Gl sparticles will be discovered and some of their properties measured. Hopefully, the 
lightest Higgs scalar h could be discovered, and some of the squarks and the gluino could be copiously produced 
(if not too heavy) at the LHC due to their strong interactions. In addition some of the neutralinos, including the 
lightest supersymmetric sparticle (LSP), could be identified and have their masses extracted indirectly from a detailed 
study of squark and gluino cascade decaysJi|i01. Beyond that, the discovery and measurement of the full set of 
MSSM sparticles may be very model dependent, and anyway challenging in many scenarios at the LHC. Various 
analyses have been conducted (see e.g. |5j,l6]) to determine the basic MSSM parameter space from the above assumed 
experimental measurements. A widely illustrated strategy, in a so-called "top-down" approach, is to start from a 
given supersymmetry-breaking model at a very high grand unification scale; predict for given input parameter values 
the superpartner spectrum at experimentally accessible energy scales; then fit this spectrum (together with possibly 
other observables like cross-sections etc) to the data in order to extract constraints on the basic model parameters. 
Constraints from past and present collider and non-collider data, with consequent prospects for the LHC and future 
linear collider (ILC), have been analyzed typically from systematic scanning of MSSM parameter space ItHs!] (though 
mostly in the constrained minimal supergravity (mSUGRA) |0] case). In addition, more elaborated \ 2 fitting procedure 
(or some generalizations! 10, l"l]|) have been also used in many such studies, together with Monte-Carlo or other 



process simulation tools Il2lll3l.il 411 . as well as other specific codes for parameter determination 1 1 5l I 1 611 . On general 



grounds, fitting and minimization procedures are efficient when the number of independent measurements is (much) 
greater than the number of fitted parameters of the underlying model, and provided that data are reasonably accurate. 
But clearly the minimization becomes less controllable 1 in a general MSSM with more than ~ 20 relevant basic 
parameters (even when neglecting flavor mixing in the sfermion sector). Alternatively, so-called inverse or bottom 
up reconstruction approaches are often motivated!^, 18, 19, 20 



21 



22. 



2311 . Also, a growing number of analysis 
for the LHC or the ILC appeared recently, attempting to go beyond conventional top-down fitting techniques||23J or 
supplementing these with more elaborated frequentist or Bayesian methods, with Markov chain Monte-Carlo (MCMC) 
techniques [24] in particular ifl \v\ Eill . Yet it has been stressed (for instance in refs. |23|) that the mapping from LHC 
data to the underlying basic MSSM parameters may be far from unique. However, most works still rely essentially on 
simulation tools fed with top-down MSSM Lagrangian-to-spectrum relations, while to our knowledge reconstruction 



scenarios based on explicitly inverted relations (see e.g. 1118 



,12111 ) appear not so widely explored in the literature. 
Ilql often considered rather optimistic LHC or ILC 



Moreover, many studies on MSSM parameter reconstruction! 5, 6 
scenarios, in the sense that the results presented were obtained by assuming that most, if not all, MSSM sparticles 
masses and other relevant observables are measured with the best expected accuracy. At the same time it is often 
assumed that the more constrained mSUGRA model Q (with four continuous plus one discrete parameters) is to be 
determined. While such studies are certainly very useful guidelines for LHC and ILC analysis, these assumptions may 
be considered quite optimistic for the supersymmetry discovery prospects in general, especially at the LHC. It is thus 
worth to develop alternative (or rather complementary) strategies to foresee more pessimistic scenarios, still trying 
to extract as much as possible informations on the nature of the underlying supersymmetry-breaking model in case 
where only a handful of the predicted sparticles would be identified. 

To start with, in this paper we explore specific bottom-up reconstructions, which are more restricted and certainly 
far from being fully realistic as concerns data simulations, but that we expect to be useful and complementary to 
the more standard simulation tools. Our approach is based essentially on analytical inverse relations between the 



1 See however ref. fnll for a recent elaborated treatment of many-parameter cases. 



measured masses and basic parameters. This "inverse mapping" for the MSSM spectrum has been investigated to 1 
some extent in the past years lll8l Il9ll but mainly at the tree-level approximation, and moreover much often in the 
context of the ILC data essentially. It is generally expected that simple analytic expressions between observables and 
parameters are more transparent or insightful than purely numerical results, providing e.g. explicit correlations among 
parameters. Particularly in the MSSM, even at tree-level, this connexion is already quite involved so that it is difficult 
to grasp with a good intuition the sensitivity of the different observables to MSSM parameters, unless having spend 
much time in doing fits and related calculations. But more concretely than a useful insight, we also hope that such 
an approach could suggest new strategies for reconstruction of parameters, as will be illustrated here. For example, 
by exploiting well-known relations between the (first two generation) squark and slep ton soft mass parameters and 
physical masses, including the renormalization group evolution (RGE) dependence 11261 12711 . we construct appropriate 
combinations of observables in this squark/slepton sector which appear to provide interesting and almost model- 
independent constraints. 

Deriving analytic inversion relations in the MSSM may appear at first a rather academic exercise, quite remote from 
the complexity of the actual experimental situation especially at the LHC. This is because such inverted mapping 
remains relatively simple only if restricted to the tree-level approximation. But it becomes a priori inextricable if 
including radiative corrections, which are certainly necessary at the accuracy level expected for realistic LHC and ILC 
data analysis. More precisely through loop contributions almost all sparticle masses have a cumbersome dependence 
on almost all MSSM parameters. Still, we will see how radiative corrections can be incorporated into our framework 
rather simply, essentially by (numerical) iterative procedure in a manageable way, in reasonable but often realistic 
approximations. We emphasize that this procedure is very similar to the way in which radiative corrections are included 



in more conventional top-down MSSM spectrum calculations l28[ l29l |30L 13 ill , and it allows to keep most advantages 
of the bottom-up approach. 

Even if one can incorporate a fair amount of presently know radiative corrections into this framework, we stress 
that our motivation is not to compete with the state of the art in present analysis of MSSM constraints at LHC, merely 
by replacing elaborated simulations tools with a bunch of rather simple analytic relations (and simple combinations 
of data uncertainties as we will see). Accordingly our analysis at this stage is still essentially a theoretical exercise, 
not incorporating important ingredients of the complexity of LHC measurements (such as detailed event selections, 
detector properties etc) that are ultimately necessary and that pave the non-trivial steps in going from basic LHC 
data to sparticle mass measurements. Yet our aim is to consider as much as possible realistic and minimal LHC 
sparticle identifications, using a limited set of sparticle mass measurements. We then gradually consider different 
SUSY-discovery scenarios, going from minimal input assumptions to more optimistic ones, defining corresponding 
algorithms with definite input/output parameters. This step-by-step analysis may turn out to be closer to the actual 
experimental situation, in which one will certainly not identify all sparticles at once, even for the most optimistic 
expectations. However our approach is only one step in the vastly more ambitious program of so-called "blind" 
analysis of LHC data: in particular when we consider a general MSSM case (i.e. departing from e.g. mSUGRA 
universality relations), we nevertheless assume a spectrum pattern still allowing gluino/squark (long) cascade decays 
(i.e. with some neutralinos sufficiently heavy to decay in the LSP plus sfermions, but sufficiently light to be decay 
product of heavier squarks/gluinos). This pattern may admittedly be considered a not so general scenario within the 
MSSM. 

Though our analysis essentially concentrates on sparticles expected to be accessible at the LHC, it will appear that 
some of the reconstruction algorithms used here could apply more or less directly to ILC measurements, upon 
appropriate changes in data accuracies. We thus occasionally make some comments on ILC prospects, but refrain to 
enter a detailed ILC analysis which is beyond the scope of the present paper, since the other expected sparticle mass 
measurements at the ILC would need rather different algorithms (though quite similar in spirit). 



The paper is organized as follows: in section 2 we briefly define and review a plausible set of sparticle mass mea- 



surements at the LHC, with accuracy on which is based our analysis. We consider different levels of assumptions on" 
the nature and number of identified sparticles, defining several scenarios. We also gradually introduce universality 
assumptions for the soft-SUSY breaking parameters of the different sectors. In sections 3-6 we expand results of 
analytic inversion analysis, with some of these already presented in ref. 11811 . for different parameter sectors of the 
MSSM, recasting results in the context of gluino/squark cascade decay mass measurements at the LHC, and incor- 
porating radiative corrections. We consider separately four different sectors: gaugino/Higgsino parameters (section 
3); squarks and sleptons (first and second generation) (section 4); third generation squarks (section 5); and finally the 
Higgs parameter sector in section 6. These distinctions are quite natural when considering both the interdependence 
between observables and parameters and the experimental signatures which are expected at LHC from a given sector. 
We will delineate which relations and results are valid in a general (unconstrained) or a more constrained MSSM (with 
universality relations at the GUT scale). We also compare in some detail at different stages our reconstruction results 
with more standard top-down fitting procedure using MINUIT x 2 minimization ! 1 ill , with data and fitted parameters 
directly set by the above step-by-step scenarios, rather than by performing "all at once" global fits. Conclusion and 
outlook are given in section 7. 

Finally we develop in Appendix A the explicit inverse solutions in the gaugino/Higgsino sector for different in- 
put/output assumptions and related issues, and in Appendix B the properties of the bottom-up renormalization group 
evolution, a necessary ingredient in this approach, implemented as an option of the SuSpect code 13 ill . Important 
features, like the error propagation from low to high energy parameters that is implied by RGE, are illustrated there. 



II. BOTTOM-UP STRATEGY FROM PLAUSIBLE LHC MEASUREMENTS 

At the LHC, the dominant production of pairs of gluinos or squarks (or gluinos associated with squarks) is expected 
due to their strong interaction. The corresponding cross-sections are large for moderate masses but decrease rapidly 
for large gluino and/or squark masses. Discovery prospects for gluinos and squarks with masses up to a few TeVs 
are reported |il EE depending on the luminosity (and depending of course on the details of the supersymmetric 
models and spectra). 

A. Mass measurements from gluino/squark cascade decays 
From detailed studies of gluino/squark cascade decay products at the LHC, the masses of the sparticles involved can 



be determined with a quite good accuracy (a few percent) using the so-called kinematic endpoints method[3, 4]. For 
a typical mSUGRA benchmark point like SPSla [34J], which has been intensively studied in simulations, the masses 
of the sparticle involved in the gluino and squark decays are obtained from analysis of exclusive chain of (two-body) 
cascade decays, typically 01 : 

9 ~> qL{b)q{b) -> x° Ifl ~> W<7/<7 -> Xi Ifkfi ■ (2-1) 

Note however that we will subsequently use these data as a blind input, with the aim to go beyond the SPS 1 a benchmark 
(or even beyond a mSUGRA model) as concerns the basic MSSM parameter reconstruction. 

Actually the four masses of gx, cr and Xi 2 (designated in what follows as , rrijj ) can be determined from the 
cascade decay starting from a qx- The gluino mass nig can then be determined from the decay to b (see B. Gjelsten 
et al p. 213 in 0]). Alternatively in ref. the gluino mass as well as the four other masses are determined from 
the full cascade Eq. d2.lt . The sparticle mass determinations and accuracies assumed here are based on the results of 
ref. together with ref. 0, summarized in Tables HI and Hill These accuracies in Table HTTl may be subject to some 
adjustments or updates due to eventually more refined analysis, and should be considered here as illustrative, without 



drastically changin g our p rocedure and results. (For instance very recently even better prospects on mass accuracie^ 
have been reported [35l|36[] by exploiting correlated decays of two such cascades.) 

Among the different selection criteria, an important characteristic one is to look for two isolated, opposite sign, 
same flavour e, /i leptonsfl Etl- We will not be involved here with a concrete analysis of these cascade events, and 
rather use directly the expected mass measurements extracted from such studies. We refer to these references for 
more details and shall only briefly mention some of their main features. For typical benchmark points, like SPSla or 
other cases[34], more sparticles than those present in Eq. ( 12. Il l are in principle accessible from independent processes. 
(Indeed the slepton Ir can also be measured independently of the cascade ( 12. U from slepton pair production). Some 
of the other sparticles may be more difficult to identify, due e.g. to the fact that neutralinos, and charginos decay 
predominantly into f and r, experimentally more challenging to detect than a dilepton signal typically ||5|]. (These 
effects are more pronounced for large tan (3 due to a larger mixing). Though gluinos decay predominantly in 
squark, and xi predominantly in Ir slepton, the mass measurements of qR (and l£) are also possible^, but may be 
less favored by the small g — > q^q branching ratio (B.R.) with respect to other channels, though final statistics can be 
sufficient fl]. Moreover qR decays directly into the lightest neutralino xi> since xi is mainly Bino for SPSla (and this 
is more or less so in most mSUGRA cases as well, similarly the next-to-lightest neutralino %2 is essentially Wino). 
Thus, qi-/ being a SU(2) singlet, it decays into the corresponding qR quark together with xi with a B.R. of almost 
100%. 



TABLE I: different plausible gradually optimistic assumptions on the amount of sparticle mass measurements at the LHC, from 
gluino cascade and other decays, defining our input for various reconstruction scenarios. 



Input scenarios 


mass 


decay or process 


(+theory assumptions) 






(minimal): 


m- g , 


g cascade decay 


Si (MSSM), 






& (universality) 


m N 2 - 




Si, 


niq L , 




S'4 (universality) 


mr 




S3 = Si plus: 


m iV 4 


qh —> xi + • ■ cascade 


S B . 




g cascade decay 


S5 (universality) 






S 6 = S 2 + S'i + S's plus: 


m h 


h — > 77 (mainly) 



More generally the nature of sparticles involved in such cascades or other considered processes strongly depends 
on MSSM parameters, i.e. the specific masses, branching ratios and other observables, and also on some properties 
inherent to the MSSM structure. Thus at present it is hard to guess which process may be actually favoured at the 
LHC. Indeed the parameter space where decay chains such as in Eq. d2.11 i can occur may be considered already 
quite specific, as it requires (non-LSP) neutralinos heavier than the first two generation sleptons, but light enough to 
be produced by gluinos and squarks. Accordingly we insist that considering in this work only the sparticle masses 
accessible from the decay d2.lt (plus the lightest Higgs) is not a strong prejudice against the possibilities of other 
processes and extra sparticle identification. As motivated in introduction, it is simply to consider what this approach 
can do from a well-defined "minimal" input set, and indeed most of our inversion algorithms could be easily extended 
if more (or different) sparticles will be available. 

Another important feature of the decay in Eq. ( 12.1b is that there is no way to distinguish the different squarks q 
from each others: this is not so much a property of this specific decay but rather due to the fact that there is no realistic 



mean at present of tagging light quark charge and/or flavor (moreover they all have almost indistinguishable B.R.^ 
Accordingly the first squark entering the decay chain, resulting from the decaying gluino, can be either u, d,c, s or 
b\ , 62 (in general it could also be t\ , but this is not kinematically allowed for the SPS 1 a input parameters Jj, |5[]). 
One can identify the b to some extent: the decay of a gluino into b\ is dominant over the 62 one due to the smaller 
mass, and the b\ decay leads to a b-quark that can be tagged. In addition, one may be able to extract a signal even for 
62 (i.e. distinguish it from 61), but with less statistics (correspondingly with a larger mass error), and only for the large 
luminosity prospect of 300 fb -1 !^. We will thus consider in addition to our minimal input scenario a next scenario 
where either 61 alone or both 61 and 62 masses can be extracted. 



Finally, on top of the sparticle masses measured via the gluino cascade, we will consider in section 6 what additional 
constraints are obtained within our approach if the lightest Higgs mass is assumed to be measured via its 77 decay 
modes Jil, which is mainly responsible of the expected accuracy as quoted in Table HIH 



B. Outline of bottom-up reconstruction algorithms 



According to the previous experimental possibilities, we define in Table Q] successive scenarios to be studied and 
differing on the amount of sparticle masses measured at the LHC, from S\ to Se- Scenarios S1-S3 may be considered 
to range from a minimal input scenario to gradually more optimistic ones, while some scenarios differ by model 
assumptions (general MSSM, or with universality relations in the sfermion and/or gaugino sectors typically). 

In our study we shall first generate "data" with central values, e.g. for the SPS la benchmark point, by running the 
code SuSpectj3l|] for the (constrained MSSM) input: 

m = -A) = 100 GeV, m 1/2 = 250GeV, tan/3(m z ) = 10, fi > 0. (2.2) 

The resulting spectrum in TableHIlis calculated from the latest version 2.41 of SuSpect. We used two available options 
on RGE and sparticle mass radiative corrections, both for illustration and the need of our analysis, as will be developed 
later on. Note that we use everywhere a value of the top mass m t0 p = 175 GeV rather than the latest experimental top 
mass values: mt op = 172.6 ±1.4 GeV 13711 . in order to be more consistent with the analysis performed in ref. 

Hi. 



We assume that this shift in the central value of the top mass should not affect qualitatively our analysis (although 
what could be important is the impact of the top mass uncertainties). 

We then use the sparticle masses contributing to the gluino cascade as blind input, within different reconstruction 
scenarios, without necessarily assuming a constrained MSSM with universality relations. The aim is to examine what 
can be reconstructed under different gradually constrained assumptions on the MSSM parameters. We illustrate the 
determination uncertainties from the SPS la sparticle mass error as reference, since it is one of the most simulated 
benchmark in the literature. 

Apart from distinguishing different scenarios as indicated in Table |TJ most of our study is based on specific 
bottom-up algorithms depending on the assumed input sparticle masses and output basic parameters. In defining 
these algorithms it is convenient to consider separately and gradually the three different sectors of gaugino/Higgsino, 
squark/slepton, and Higgs sector respectively (distinguishing also the third from the first two generations in the squark 
sector, since those necessitate different treatments due to the mixing in the third generation). We will also carefully 
distinguish different scenarios depending on the amount of theoretical assumptions, eventually reducing the number 
of basic MSSM parameters, like universality of gaugino and/or scalar mass terms typically. Our different algorithms 
obviously depend on specific assumptions, since input and output parameters may be completely different depending 
on these. We describe in detail in the next sections these particular algorithms depending on the parameter sectors and 
theoretical assumptions considered. The starting point is always the use of tree-level relations giving some specific 
Lagrangian parameters in terms of appropriate input sparticle masses. For example in the gaugino/higgsino sector, one 
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TABLE II: Soft and other basic parameters, plus sparticle pole masses (relevant to our study) for SPSla input (with m top = 175 
GeV), calculated with the latest SuSpect ver 2.41, for two illustrative optional choices: 1) full two-loop in RGE and full radiative 
corrections to sparticle masses (second and fifth columns); 2) one-loop RGE, no radiative corrections to squarks, gluino, neutralinos, 
charginos masses, simple approximation for mj, radiative corrections (third and sixth columns). 



basic par. 


2-loop RGE 


1-loop RGE 


relevant 


2-loop RGE 


1-loop RGE 




+full R.C. 


+approx. R.C. 


pole masses 


+full R.C. 


+approx. R.C. 


Qewsb 


465.5 


468.2 








Mi 


101.5 


108.8 


m jV! 


97.2 


105.1 


M 2 


191.6 


208.9 


m N 2 


180.8 


189.9 


M 3 


586.6 


603.8 


rrig 


606.1 


603.8 


A 1 


356.9 


340.6 


IV 4 


381.8 


369.6 


tan (3 


9.74 


9.75 








11 d 


(179.9) 2 


(187.3) 2 


m h 


110.85 


111.28 


m 2 Hu 


-(358.1) 2 


-(341. 7) 2 








m eL 


195.5 


201.5 








m TL 


194.7 


200.6 








m eR 


136 


138.6 


mg 2 


142.8 


145.4 


m TR 


133.5 


136.2 








™>Q L 1,2 


545.8 


554.1 




562.3 


551.6 


m Q L 3 


497 


502.9 


m ' bl 


516.2 


502.1 


m UR 


527.8 


531.6 








m tR 


421.5 


421.6 








m dR 


525.7 


528.7 








m bR 


522.4 


525.4 




546.3 


530.1 


-A t 


494.5 


501.0 








-A b 


795.2 


791.3 








—At 


251.7 


255.0 










677.3 


686.6 








-Aa 


859.4 


857.2 








-A e 


253.4 


256.7 









of the inverted relations we shall consider has the form 

f(mfi 1 ,mff 2 ;M 1 ,M 2 ) + Af rad . corr . = 0, (2.3) 

where f gives the output parameters Mi , M2 in terms of two neutralino mass input: rrift , m# > or other such relations 
for different input/output choices. Whenever possible, the relations defining / in Eq. (12.3b are entirely analytical 
and often giving a linear (unique) or at most quadratic solution (with eventually corresponding twofold solutions). In 
addition some input/output parameter choices need extra numerical calculations, typically iterations. These are needed 
anyway to take into account the radiative corrections, symbolized by the term Af ra d. C orr. in Eq. ( 12.31 l. which generally 
depend on extra MSSM parameters or masses. As already mentioned, it is clear that such approach cannot be very 
realistic if not including at least some part of these radiative corrections, as is discussed in next sub-section. 

Once having reconstructed from a relation like ( 12.31 ) the relevant MSSM parameters at the "physical" scale (gen- 
erally identified as the electroweak symmetry breaking (EWSB) scale), another important step in this bottom-up 
approach is the possibility of evolving these parameters consistently from low to high (GUT) scale, with implications 
concerning the propagation of parameter uncertainties from low to high scales. Such bottom-up RGE evolution of soft 



parameters had been considered in the past il 18L 12011 (see 12011 notably for mass measurement error propagation), but 



meanwhile many refinements e.g. on radiative corrections have been included in public MSSM codes. Accordingly we 



have implemented an up-to-date option in the code SuSpect to perform this bottom-up RGE, which is used at differerS 
stages in our analysis and illustrated in more detail in Appendix B. 



C. Including radiative corrections in bottom-up reconstruction 



X; Yj (GUT) Xj Yj (GUT) 




FIG. 1: Schematic algorithms of the top-down (left) and bottom-up (right) mappings illustrating their similarities in practical cal- 
culations. Xi and Yj are a set of (soft or SUSY) running parameters, RG-evolved between GUT and EWSB scales (steps (1) or 
(3) respectively). Xi and Yj may mix, and diagonalization gives (running) mass eigenvalues mx , my (step (2)). Then, additional 
radiative corrections linking running to pole masses m l pol& are added (resp. subtracted) and may depend on extra unknown param- 
eters or masses Zk,rrik- Iterations are needed at several steps in each approaches, as well as specific assumptions on some a priori 
unknown parameters. 



We explain here on general grounds how we incorporate radiative corrections into our algorithms, with specific 
mass/parameter relations to be given later, once having defined algorithms for the different sectors more precisely. 
Prior to the bottom-up RGE comes the question of incorporating radiative corrections linking the running parameters 
to the physical (pole) masses, as generically indicated by the terms A/ mc j. corr , in Eq. ( 12.31 ). Clearly, incorporating 
the full radiative corrections to all sparticle pole masses irremediably spoils such simple analytic inversions, since the 
complete radiative corrections would introduce, already at one-loop level, highly non-linear dependence upon (almost) 
all parameters of the MSSM. However, upon assuming a certain level of (reasonably good) approximations for these 
radiative corrections, it turns out to be relatively easy to incorporate these at a realistic level. This is especially the 
case for the sparticle masses entering the relevant cascade decays: typically the (first two generation) squarks receive 
radiative corrections that are largely dominated 138|] (at one-loop) by gluon/squark and gluino/quarks QCD corrections, 
involving precisely the same sparticles entering the cascade. Other corrections, like the electroweak ones, are fairly 
negligible ]38|] in comparison. Similarly, corrections to the gluino mass are essentially dominated by gluon/squark 
QCD corrections. Consequently it is rather simple to "subtract out" those corrections, starting from the experimentally 
measured pole masses, and applying next tree-level algorithms to the running masses. It eventually needs to apply 
this procedure in several steps using numerical iterations. Although this procedure may appear rather involved, we 
emphasize that it is very similar to the manipulations which are p erformed in a standard top-down approach, where 
iterations are anyway necessary in spectrum calculations 31] once including radiative corrections. The 
similarities between standard top-down and our bottom-up practical calculations are illustrated sketchily in Fig.[T] Xi 
and Yj represent some of the relevant (soft or SUSY) running parameters, to be RG-evolved between GUT and EWSB 
scales (step (1) in the top-down approach or (3) in the bottom-up approach respectively). Then these parameters can 
have mixing so that diagonalization in step (2) gives (running) mass eigenvalues mx > my in the top-down case. In 
the bottom-up case, rather than performing a brute force inverse diagonalization at step (2), it is more convenient to 



use appropriate relations 1 1811 such that the required output MSSM parameters are in one-to-one relations with the few 1 
accessible input masses. (We will see a specific example of such relations for Eq. (12.3b in the neutralino sector in 
section 3 below). Next, the necessary radiative corrections which link the running to the pole masses mJ ole are added 
(respectively subtracted in the bottom-up case). These may depend on extra unknown parameters or masses Zjt, 
(in which case definite assumptions on these unknown parameters are needed). Iterations are performed in each 
approaches for the radiative correction steps, since these depend on final sparticle masses. (Iterations are also needed 
for the RGE since the EWSB scale and other relevant parameters depend on the sparticle spectrum). In practice the 
kind of subtractions and other numerical manipulations that are needed specifically here are made easier by a number 
of possibilities included in the latest version of the SuSpect lBlll code 2 . Concerning the neutralino masses, radiative 
corrections are know to be reasonably small, and moreover to a very good approximation one can incorporate the 
leading ones in the form of tree-level deviations on the parameter fx, Mi, M2, allowing again subtractions and iteration 
procedures when applying tree-level reconstruction algorithms. Moreover, in some cases we also incorporate the extra 
unknown radiative corrections by assuming typically universality relations within the loop-level calculations. This 
may induce a little bias, but we consider (and have explicitly checked for the SPSla case) this to be a reasonably good 
approximation, even when considering a general MSSM reconstruction. 

Concerning the Higgs sector, radiative corrections to the light Higgs mass rrih and pseudoscalar mass rtiA are known 
to be of primary importance. But, as is well known, the leading contributions essentially come from the stop sector, 
and more generally there exist approximations [|39l that are excellent to 1-2 GeV level 3 , i.e. of the order of higher 
order uncertainties j4oll . 

On general grounds, even within the present state of the art, the known radiative corrections to sparticle and Higgs 
mass still suffer from uncertainties due to unknown higher orders. Moreover at the LHC experimental errors are 
generally larger than the latter theoretical errors (except for the lightest Higgs mass rrih). These features evidently 
affect our analysis, but in the same way as any other more standard top-down approach to the reconstruction of MSSM 
parameters at the LHC. Clearly the real limitation in incorporating radiative corrections does not come from the 
eventual complexity of incorporating these numerically within a particular procedure, but rather on the uncertainties 
resulting from unknown sparticle masses contributing at the loop level to a given observable. In some cases where the 
latter uncertainties may particularly affect our results, we take these into account as theoretical uncertainties, as will 
be specified. Overall we consider that our treatment of radiative corrections as described here should be sufficient for 
our rather limited purpose. 



D. Treatment of mass uncertainties and interpretation 



For the different scenarios considered we will illustrate the expected accuracy on the reconstructed parameters for 
given mass measurement accuracies assumed according to Table Hill (eventually considering also theoretical errors). 
To delineate this error propagation we have performed various scanning over the input mass values within errors, 
or over relevant MSSM parameters, either with uniformly distributed random numbers, or alternatively also using 
random numbers with a Gaussian distribution (in which case we can define confidence level intervals) 4 . The sparticle 
mass errors as quoted in Table Hill are, however, known to be not purely statistical: there is a large part which comes 



2 like e.g. the option to "switch off' gradually some of the radiative corrections to the sparticle masses, as illustrated in Table ITTl 

3 The latter approximations are also incorporated as one option in the SuSpect code, alternatively to the full one-loop, or full one-loop plus leading 
two-loop calculations options. 

4 There are a few cases in our analysis where uniform "flat prior" distributions may give misleading "density" regions for the resulting constraints 
on some of the parameters (typically for tan (5, see sec. 4.2 for illustration and discussion). In such cases we made obvious changes for more 
appr opriate non-flat distributions, but have not made any attempt to define much refined priors in a Bayesian approach such as is done notably in 



TABLE III: Experimental accuracies on mass determinations from LHC gluino cascade and other decays (taken from refs.[2|] and 
@]), corresponding to the different input mass scenarios defined in Table[l] 



mass 


expected LHC 
accuracy (GeV) 


decay or process 


m- g 


7.2 


g cascade decay 




3.7 




m N 2 


3.6 






3.7 




In 


6.0 




m N A 


5.1 


Ql — > X4 + ■• cascade 


m - bl 


7.5 


g cascade decay 


m h 


7.9 




m h 


0.25 (exp)-2 (th) 


h — > 77 (mainly) 



from the systematic errors on jet resolutionyl |4J], and moreover these errors are also strongly correlated. We stress 
however that a more involved treatment of uncertainties, properly combining the statistic and systematic ones, taking 
into account correlations etc, appears quite non-trivial[17J and is beyond the scope of the present paper. One could 
in principle make substantial improvement in the final determination of parameters by using directly the endpoint 
measurements of the gluino cascade rather than the naive mass errors obtained from the latter. Consequently, 
one should keep in mind that the interpretation of the various domains and contours in parameter space that we shall 
obtain are lacking a very precise statistical significance. (We plan to perform a more refined statistical analysis in 
the future 114 ill ). Despite these limitations, we will illustrate detailed comparisons for most considered scenarios of 
bottom-up determination results with those obtained from more standard statistical treatment with % 2 minimization in 
a top-down approach using MINUITfl 111. 



III. GAUGINO/HIGGSINO PARAMETER DETERMINATION FROM GLUINO CASCADE 

We start by recalling some analytic inversion algorithms at the tree-level, adapted to the LHC input scenarios (i.e. 
corresponding to the extractable masses in the gluino decay chain as discussed above). Note that our algorithms may 
be valid more generally, e.g. at the ILC, provided that the same input masses would be available. 

A. Gaugino/Higgsino parameter: general case inversion 



We thus consider the parameters relevant to the gaugino/Higgsino sector, starting from the neutralino mass matrix: 



M N = 



Mi —nizsw cos [3 mzsw sin/3 \ 

Mi mzcw cos (3 —mzCw sin /3 

— mzsw cos/3 mzcw cos /3 — fi 

\ mzsw sin (3 —mzCw sin /3 — // ) 



(3.1) 



Rather than performing an involved inverse diagonalization, which would moreover need to know all the four neu- 
tralino masses, it is much more convenient to use appropriate relations among parameters involving fewer input 



masses. The four invariants (under diagonalization transformation): 



11 



TrM N , - Tim 



(TrMjy) 3 _ 7YM IY(Mj» ) + Tr(M|) ^ ^ 
6 2 3 



provide a system of equations] 18] which can be used in different ways depending on the choice of input and output 



parameters. Two equations are actually expressing necessary and sufficient conditions for the existence of solutions to 



this system (see also Appendix B of ref. [18] for more details) 



P 2 + {[? + m%- MiM 2 + (M, + M 2 )Sij - S 2 J )P lJ + ^m%{c 2 v M 1 + s 2 w M 2 ) sin 2/3 - /i 2 MiM 2 = (3.3) 
and 

(M 1 + M 2 - S l3 )P 2 + {n 2 {M x + M 2 ) + mKc^M, + s 2 w M 2 - M sin2/3))P y 

+ f x(m%(c 2 v M 1 + s 2 w M 2 ) sin 2/3 - nM 1 M 2 )S ij = (3.4) 

where we define for short SV, = m^.+m^., Pjj = m^.m^, where i, j = 1, ...4 5 , and s\y = sin^iy, cw = cos @w- 
Note that Eqs. ( 13.31 ), ( 13.41 ) involve only two neutralino masses, which corresponds to our minimal input assumptions 
in Table U] These are originally tree-level relations but, as explained in sub-section 2.3, in our analysis we shall 
incorporate as much as possible of realistic radiative corrections. To begin, the values of and mz in expressions 
(13 . 3 b , <f3~4b are understood as the properly defined DR scheme parameters: siy and mz- 

If chargino masses were known at this stage Eqs. ( 13.31 ). ( 13.41 ) would lead rather simply to a unique solution for 



Mi for given fi, M 2 and tan /3 lll8h . This had been studied in the past for chargino and neutralino mass measurement 



prospects at the ILC. Precise determinations of the chargino/neutralino parameters at the ILC, partly based on analytic 



(or semi-analytic) inverted relations in the neutralino and chargino sector, have been largely analysed in ref. [ 19]. But 
since we do not assume chargino masses to be measured in our scenarios (which appears anyway more challenging 
at LHC), and given the parameters entering the relevant gluino/squark decay, it is more appropriate to use Eqs. ( 13. 3b . 
( 13.4b differently as we examine next 6 . 



1. Scenario SI: determining Mi, Mi from m^, mpf 2 in general MSSM 

We first consider a general (unconstrained) MSSM scenario SI, assuming non-universality of gaugino masses. We 
then use Eqs. (13. 31 ). (13.4b to determine Mi and M 2 from (any) two neutralino mass input: we thus take m^, 
input extracted from the cascade decay, for given /i and tan (3 parameter input. It is straightforward after some algebra 
to work out from Eqs. d3.3l ), ( 13.4b these Mi,M 2 solutions (e.g. eliminating first Mi which depends linearly on M 2 
from one of the two equations, and obtaining a quadratic equation for M 2 ). For completeness the explicit solutions 
and related issues are worked out in some detail in Appendix A (see Eqs. ( lAlb - dA3l ). We note here that the solution 
for Mi, M 2 has actually a twofold ambiguity, being obtained from a quadratic equation e.g. for M 2 . More basically 
Eqs. ( 13.3b . ( 13. 4b as well as all other relations (13.2b only use information on mass eigenvalues, and are invariant under 
any neutralino mass permutations, e.g. <-» m^. Accordingly, without further theoretical assumptions on 

gaugino mass terms, one cannot establish the hierarchy between the two gaugino (and the Higgsino) mass parameters 



5 Note that , can be any two neutralinos, all these equations being symmetrical under any neutralino mass permutations. 

6 NB another recent analysis of the neutralino system in the LHC context of gluino/squark cascade decays has been performed in ref. |42|| , also 
partly based on semi-analytic relations, though very different from ours and not relying on exactly the same input. 



from the sole knowledge of those two neutralino masses, unless extra information on the diagonalizing matrix elemerki 
is available (which amounts to have information on some of the neutralino couplings to other particles). Thus in a 
general gaugino mass scenario there are two cases to consider, depending on the relative values of the Bino and Wino 
soft mass terms: either Mi < M<x, as in most mSUGRA scenarios, or a reverse hierarchy My, < Mi (as in the case of 
e.g. AMSB models). 

When assuming a well-defined Bino/Wino mass hierarchy, the (Mi, M2) solution is then unique 1 . Now taking 
central values of the masses 97i# , rrijj plus the reference SPSla values of fi and tan (3 we recover the correct SPSla 
values of Mi, M2 if assuming Mi < M2, or another possible solution in general MSSM with M2 < Mi as is 
examined further below (see also Appendix A for more details). More interesting than this explicit solution for 
fixed input values is to determine the expected accuracy on output parameters, given the experimental uncertainties 
on neutralino masses, and the sensitivity of Mi , M2 to the presumably limited knowledge on the two other basic 
parameters /1 and tan (3. This error propagation and other issues in the reconstruction of Mi , M2 for the SPS 1 a test 
case will be illustrated below in subsection lHIBI 



2. Scenario S2: determining fj,, tan (3 with gaugino mass universality 

In a different scenario we consider the very same basic Eqs. (13.31 1. (13. 4t but changing input/output: adding now the 
gaugino universality assumption: Mi = M2 = M3 at the GUT scale, we first determine Mi, M2 from M3, at the 
EWSB scale. (This does not necessarily imply a mSUGRA model, since non-universal relations could still hold for 
all other MSSM parameters apart gaugino mass terms. At this stage one could also start from any other well-defined 
relation between the M^'s at some given scale, like is the case for AMSB and GMSB models). As a consequence of 
the related RGE structure of gaugino masses and gauge couplings at one-loop level, the relation in the universality 
case reads: 

Mi Ah _ M3 n - 

9 9 9 

91 9 Z 2 33 

(where are the properly normalized gauge couplings) to be valid at any scale. Then, Eqs. ( 13-31 ), ( 13.41 ) are now used 
to determine fj, and tan/? for (universal) Mi, M2 input, as a linear system for sin 2(3 and /1 2 . It is simple after some 
algebra to work out those explicit solutions, e.g. first eliminating sin 2(3 to get an expression for /1 2 that only depends 
on Mi, M2, and the two input neutralino masses. For completeness explicit solutions are given in Appendix A (see 
Eqs. ( lA7b . (1A9I >). Our conventions are the usual ones such that < (3 < ir/2, so that tan/3 > (and real). This is 
not a restriction on parameter space, since an eventual phase of tan (3 can be absorbed by a consistent redefinition of 



the Higgs doublet fields 12211 . The sign of fi, however, is not determined by these equations, so we have to consider the 
two possible solutions for [i > and /1 < a priori. As previously, as a cross-check we can plug in these expressions 
the central SPSla values for the masses M§, Mjvi, Mm 2 as obtained e.g from SuSpect, obtaining the correct values 
of tan (3 and fi. In the numerical applications for SPSla reconstruction, illustrated in subsection lHIBI below. we shall 
thus consider both \i > or \i < case (examining whether the latter may be eventually eliminated when taking into 
account input mass accuracies.) 



3. Scenario S3: three neutralino mass input (with and without gaugino universality) 

What could be more constraining is the (more optimistic) scenario where three neutralino masses could be deter- 
mined at the LHC, involving another squark decay measurement (independent from the first gluino cascade decay)!^, 



7 We assume M2 > in addition, which one always has the freedom to choose in MSSM l23l . 



according to the input S3 in Table Q] above. In this case one can use very simply an extra relation originating frcM 
Eqs. ( 13.2b to get a determination of either (x or tan/3. More precisely from the trace of the matrix (13. lb and the second 
invariant in Eqs. ( 13.2b , one obtains a simple expression for /j 2 : 

= M X M 2 - m| - (F 124 + SiM{Mi + M 2 - S 124 )) (3.6) 

where S124 = + mjy- +111^^ and P124 = m^jn^ + rnfj mfj + . Eq. (13.6b can be first used 

in the non-universal scenario S\ above, thus determining Mi, M2 and from three neutralino masses (plus tan/3) 
input. (Alternatively one may also solve this system for tan ft instead of \i, but since all expressions only depend 
on sin 2/3, it becomes rapidly insensitive for large enough tan ft. Accordingly we anticipate without calculations that 
it is unlikely to get any interesting tan ft (upper) bounds given the input mass LHC accuracies, irrespectively of the 
amount of neutralino masses measured.) Solving Eq. ( 13.6b together with Eqs. (13.3b . (13.4b gives in fact a high (sixth) 
order polynomial equation for Mi, M2, /x which thus cannot be solved fully analytically. It is however easy to solve 
iteratively using e.g. Eq. ( 13. 6b on the solutions ( IA2b , ( lAlb (upon having chosen a definite Mi, M2 hierarchy). This 
iterative solution converges very quickly (see Appendix A for more details). 

When applied to the reconstruction for SPSla test example, with corresponding input mass error propagation, this 
will result in a much more precise determination of fj,, as will be illustrated in subsection lHIDI Note however that the 
sign of fj, remains undetermined from this additional information. 

B. Reconstructing Mi, M2 in MSSM without universality assumptions: SPSla test case 

We now apply the general solutions obtained for Mi, M2 in the non-universal gaugino mass scenario SI (as de- 
scribed in subsection IIII A 1 1 and detailed in Appendix A), to the actual reconstruction of those parameters for the 
SPSla test case, taking into account input mass error propagation 8 . This is shown in Fig. [2] where domains in the Mi, 
M2 plane are obtained for accuracies on the two neutralino masses taken from Table iHll resulting from a scan with 
uniformly distributed random numbers. This illustrates in particular the twofold ambiguity in reconstructing Mi, M2 
from the sole knowledge of two neutralino masses, as discussed in subsection IIII A II We also consider different as- 
sumptions on the fj, or tan ft range of variation, anticipating the difficulty in determining tan ft at the LHC solely from 
this cascade decay information, as will be confirmed more quantitatively in next sections. In practice we vary widely 
tan/3, 1 < tan/3 < 50. 

We thus illustrate the cases where both // and tan ft would be largely undetermined, and how the Mi, M2 determi- 
nation is improving if a more precise determination of fj, can be available (anticipating the better accuracies that may 
be obtained from more theoretical assumptions, or other LHC processes, or alternatively if supplementing our analysis 
with ILC determination of parameters). One observes from Fig.|2]that for largely unknown fj, and tan ft the constraint 
obtained on Ml from using solely the two neutralino mass input is fairly reasonable: 80 GeV < Mi < 120 GeV (for 
the All < M2 mSUGPvA-like pattern), and a similar accuracy for the Mi > M2 case. In contrast A/2 appears more 
poorly constrained in both cases. Moreover note that only the region M2 < 400 GeV is shown on this plot, while ac- 
tually there are a few isolated points obtained from the scan with higher M2 values, for reasons to be explained below. 
We have checked that using a regular grid scan, instead of random numbers, does not significantly alter the contours 
in Figs.[2][3]or similar other figures as will be presented below. One should indeed be careful in the interpretation of 
the density levels of various regions, since our scan was performed here with uniformly distributed random numbers. 



8 Strictly speaking one should include here the theoretical uncertainties on as well, but the latter are rather small in comparison to the 
experimental uncertainties on the other parameters. Yet, this is related to the consistent inclusion of radiative corrections, which contribute to 
5^. In fact this induces a small shift of the central values but will affect very little the variation range of the output parameters Mi, M2 here. 



Accordingly the density levels of points as appearing in Fig. [2] principally reflect that the determination of Mi, 7lif 
from Eqs. d3.3b - (l3.4b is very non-linear with respect to \i (and with respect to tan/3 to some extent), see Eqs. (IAU , 
dA2b in Appendix A, and have thus no direct meaning of statistical confidence levels. In next sections we often make 
explicit comparisons between uniform "flat prior" and Gaussian scanning of parameters: in the latter case, statistical 
confidence levels can be more properly defined (with the cautions however mentioned in sub-section 2.4, regarding 
the fact that the data used in the present work are not purely statistical anyway). In some cases the differences are 
significant and deserve a more careful analysis, as we will see. 




FIG. 2: Mi, Mi (in GeV units and at Qewsb scale) determination from two (resp. three) neutralino masses , (m^ ) in 
unconstrained MSSM for different assumptions: 1) green regions: rn^ , input, /i = /i(SPSla) ± 1 TeV, 1 < tan/3 < 50; 2) 
blue regions: nift , rrifj input, 1 < tan/3 < 50, A/i = 100 GeV; 3) red regions: m^, m^y + input, 1 < tan/3 < 50 (/j, 
is thus more constrained due to the third neutralino mass, see main text). Also shown are the mSUGRA or GMSB (resp. AMSB) 
Mi I Mi relations in orange (respectively in maroon) including experimental errors on the masses. 

The cases of moderate (blue region) and accurate (red region) determination of fj, is giving much more interesting 
constraints. The red regions are anticipating the resulting accuracy on /x, M 1; M 2 when a third neutralino can 
be measured, as will be analyzed in a next subsection IIII Dl According to Eq. ( 13.61 ) in this case /i is determined 
(independently of tan/3) and can be combined with the previous Eqs. ( 13.31 ). ( 13.4b to obtain much improved M\,M2 
determination. (Alternatively another rather good determination of fx is also obtained when the latter is not arbitrary, as 
is assumed in a general MSSM, but calculated from EWSB consistency conditions from universal Higgs and sfermion 




FIG. 3: Correlated values of A'h (left) and Mi (right) for the Mi < Mi case at Qewsb scale as functions of fi, for mg ll ,mjy 2 
input. The spreading of points is due to the variation 1 < tan/3 < 50 (in green) as well as the variation of m^, within 
accuracy (in red). 

mass terms at the GUT scale, as will be analyzed in section 6.) 9 Now, contour plots like those in Fig. [2] are not very 
informative as concerns the improvement in Mi or A/2 determination to be expected when increasing \i (or eventually 
tan/3) accuracies respectively. To trace more clearly this behaviour, we plot in Fig. [3] the equivalent of the green 
contour of Fig.[2]but in the (/i, M-x) and (/i, Mi) planes respectively (and for the M\ < M-x case). (Note that the 
values of M\ and M2 are entirely correlated since both are obtained from Eqs. ( IAU . IA2I )). The spreading of points in 
these plots is due to the variation of tan (3 and , . More precisely, what is shown in red in the effect of m ^ , 
experimental errors only, for fixed SPSla value of tan/3 = 9.74, while the additional green points correspond 
to 1 < tan/3 < 50. These plots are thus essentially the solutions Mi^/tt) from Eqs. dAU - dA2b , that would reduce 
to simple curves for fixed tan/3, m^- , m/j- . One can see the structure of solutions for (M2, /x) (and correspondingly 
(Mi,n)) with different domains, originating from the \i dependence in Eqs. iA2\ . with strong correlations. In fact 
Ma(/i) becomes arbitrarily large for two fj, values, for /j, < and fj, > (which are not exactly symmetrical, see 
Appendix A): for instance for SPSla values of tan/3, , mjy , the positive "pole" is at yu ~ 118 GeV. This explains 
the loose determination of A/2 for large variations of fi, and also explains the density levels of scanned points in Fig. [2] 
In contrast, Mi always remains finite when M2 becomes arbitrary large, according to Eq. tA.ll . This also explains the 
much better constraints on Mi in Fig. [2] irrespectively of the M2 behaviour. 

Next, one can see that both A/2 and Mi can be much better constrained, irrespectively of tan/3 values, as soon as 
the (j, determination is slightly better (such that (j, remains sufficiently far from these poles). This explains the much 
improved constraints on M2 and Mi for the blue contour in Fig. [2] (and a fortiori for the red contours where fj, is 
tightly constrained from the third neutralino mass as will be discussed more in sub-sec. UlIDI i. All these properties are 
rather simple consequences of the basic Eqs. ( 13 .3b - (l3 ,4b . and illustrate useful informations that would be very difficult 
to delineate from a more standard top-down fit of parameters. Actually the poles for specific fi values are artifacts of 
our inversion equations, but more physically it simply means that to obtain precisely the m^,m^ SPSla values for 



9 For given [i input, the measurement of a third neutralino mass could in principle resolve the twofold Mi , M2 ambiguity, since the two different 
solutions give different wijy 1 m jv 4 values. But the latter masses being essentially determined by (at least when > M\,M%), the 
differences corresponding to those two solutions are often small (e.g. only ~ 4 GeV for m ^ , for fixed SPS 1 a values of fi, which is smaller than 
the expected LHC accuracy.) So one would need to determine (and n) accurately to really disentangle the two Mi, M2 solutions. 



those particular fi values, M 2 would have to be unreasonably large. Going back to the standard top-down approach 
it also means that performing e.g. a \ 2 fit of the neutralino masses is likely to give a very flat behaviour of the % 2 
near this /j, region: more precisely, since M 2 varies widely around these \i values, no clear "best fit" M 2 value will 
be found, or with a very large error, and/or that the Xmin vame will be bad. This is fully confirmed by the results of 
a MINUIT fit: if fi is fixed to ~ 120 GeV the minimization does not give useful constraints, MINUIT finds typically 
errors like: 

Mi = 2000 ± 258 GeV, M 2 = 2000 ± 88 GeV . (3.7) 
with even many more minima and errors found for very large M 2 value. 

In contrast, fixing /j, (and tan/3) to their SPSla values and fitting only rripf , mjj , gives very good accuracy on 
Mi, M2, as will be discussed in a next sub-section below where other MINUIT fit results are given (see Table II V|>: 

Mi = 108.8 ± 5.8 GeV, M 2 = 208.9 ± 6.4 GeV . (3.8) 

Next, if gaugino mass universality at the GUT scale is assumed, as expected one obtains stronger constraints. This 
is illustrated on Fig.[2]by the (orange) band resulting from the "mSUGRA" relation: Mi ~ 0.53M2 at the low energy 
scale Q ~ Qewsb, from M\{Qgut) = M2 (Qgut) at GUT scale. The width of this band results from the error on 
M3, i.e. nig determination. We will see in next subsection how to make this study more precise when the (mSUGRA) 
gaugino mass universality is assumed. We anticipate, however, that for the given two neutralino and gluino mass ac- 
curacies, constraints on /j, will be mild, even with gaugino mass universality assumptions, while those on tan (3 almost 
absent. Indeed, one can see on Fig.[2]that the "mSUGRA" band is compatible with a part of the green region where 
tan (3 (and fi) are essentially undetermined. 

Next, since the contours in Fig. [2] are valid for arbitrary gaugino masses, it is straightforward to superpose different 
gaugino mass relations, for instance in AMSB ll^ models where the M1/M2 relation is also fixed from high scale 
boundary conditions and RG evolution, but is very different: M\(Qewsb) — %-9M2(Qewsb)- We show similarly 
on Fig. [2] this AMSB" band (in maroon) including its width originating from the M3 uncertainty. In this way one 
may possibly distinguish, depending on the accuracy, between e.g. mSUGRA/GMSB and AMSB models from the 
neutralino mass measurements. (Note however that the relation between Mi and M2 in GMSB modelsj^l is com- 
pletely indistinguishable from the mSUGRA relation at this accuracy level). More precisely one can see here how 
AMSB would be excluded if moderate (in blue) or accurate (in red) fj, measurements could be achieved (even when 
considering the second solution with Mi > M2 which has an AMSB-like hierarchy pattern). This is also a consistency 
cross-check, in the present analysis, since we started from a mSUGRA model SPSla "data". 



C. Reconstructing /i, tan j3 with gaugino mass universality for SPSla case 

Let us now consider scenario S2 as discussed above, assuming gaugino mass universality Eq. ( 13.5b to determine Mi 
and M2, and next using basic Eqs. ( 13.31 ), d3.4t to determine fi 2 and tan (3 for the SPSla test case (see explicit solutions 
Eqs. dA7t , dA8b in Appendix A). Numerically, for SPSla, gaugino mass universality at GUT scale gives at the relevant 
low energy EWSB scale approximately: 

Mi ~ 0.17M 3 , M 2 ~ 0.33M 3 , Mi ~ 0.53M 2 . (3.9) 

To determine Mi, M 2 from M3, we first extract M3 from the the gluino pole mass nig as 

irP g ole = M 3 (Q) + £§(Q;m 5 , • • •) (3.10) 

by subtracting out the leading radiative corrections Eg to the gluino mass: those are dominantly due to squarks, and 
thus largely predictable in our framework, as discussed above in sub-section 2.3. This induces a non-negligible shift, 



since for SPS 1 a the correction T, g ~ 20 Ge V, with M 3 (QewSb ) ~ 600 Ge V. 1 ' 

As mentioned in sub-sec tion llll A 2] The solutions of ( 13.31 ). ( 13.4b for tan /? and fi do not determine the sign of /i, so 
in the reconstruction with error propagation we have to consider the two possible solutions for \i > and fj, < 0. We 
first vary M g , m^, within accuracy according to Table UTTl Scanning the values with (uniformly distributed) 
random numbers, with the conditions: fi real and tan j3 > 0, is illustrated in Figs. [4}|6] A number of remarks are 
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FIG. 4: Allowed domain for corresponding to complex or real fi. 




FIG. 5: Constraints on fi > 0, tan /3 from gaugino universality and gluino plus two neutralino mass measurements in gluino 
cascade decay. 
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FIG. 6: Constraints on /i < 0, tan/3 (at Qewsb scale) from gaugino universality and gluino plus two neutralino mass measure- 
ments in gluino cascade decay. 



worth here: first, requiring \i to be real further reduces the errors on neutralino masses, since the red domain in Fig. |4] 
is to be excluded. In the most general MSSM case, /i may have a non-zero phase, but if we restrict our analysis to 
a real parameter space, this is an interesting additional constraint, resulting solely from the consistency of Eqs. d3.3t , 
d3.4t . The two other domains correspond to \i > and /i < respectively, so that the latter cannot be excluded given 
these SPSla accuracies on the two neutralino masses. We notice that if neutralino mass accuracies could be reduced 
by a factor of about 2, the \i < solution would disappear altogether (as well as the complex \x possibility). 

Next the corresponding constraint in the tan f3, fi plane are shown in Figs.[5]and|6]respectively for fx > and /i < 0. 
We observe that tan (3 is practically unconstrained, especially for large tan (3, and [i > poorly constrained, for these 
accuracies on neutralino and gluino masses. However for fj, > the two parameters appear strongly correlated, as 
shown by the contour shape: e.g. for [i > 0, large [i ~ 600 — 700 GeV is only possible for small tan j3 ~ 2. This 
correlation is not an artifact of our simple random scan, but a simple consequence of the fi and sin 2/3 dependence 
within Eqs. (13. 3t , d3.4t (see explicitly Eq. JA9b in Appendix A). Note also that the sign of /i and of sin 2/3 are partly 
correlated (see Eq. dA9K so that our convention tan (3 > impose constraints that are quite different for the [i > 
and /i < cases. More precisely we find: 

tan/3>l, 280 GeV < [i < 750 GeV (3.11) 

for fi > and 

tan/3>2.7, 210 GeV < < 370 GeV (3.12) 

for jtt < 0, and in both cases no upper limits on tan (3, which is simply understandable because only sin 2/3 appears in 
any of the relations above in Eqs. (13.2J43.4| i, so that for large tan (3 any sensitivity on tan j3 disappears. 

These results are thus consistent with what was anticipated from the previous analysis illustrated in Fig. [2] where the 
gaugino "mSUGRA" universality band is crossing all domains of the chosen tan (3, n range of variation: in particular 
the green region where tan (3 was essentially arbitrary. The fact that [i is better constrained for /i < is understandable 
from Fig. |4] where the domain corresponding to /i < is smaller than the /i > one, moreover the central value 



\n\ ~ 357 GeV is excluded on the plot Fig. [6]This is due to the partly correlated sign of fj, and sin 2(3 in this inverted 
determination (see Eq. (1A91 >), but without knowing the true SPSla value of (i we could not exclude /i < solutions 
solely from these cascade decay mass accuracies. 



D. Scenario S3: fi, tan (3 from three neutralino with or without universal gaugino masses 



Let us finally consider another (more optimistic) scenario S3 where three neutralino masses could be determined 
according to the input S3 in Table|T]above. As explained above in this case one gets from Eqs. (13 ,2b an extra relation, 
Eq. (13.6b . resulting in a determination of /i independent of tan (3, which is valid both for the general MSSM case, or 
assuming additional gaugino mass relations (like universal ones typically). For the general MSSM, the much improved 
determination of Mi, M% was illustrated by the red domains in Fig. [2] Here for completeness we illustrate in Fig. [7] 
the corresponding domains in the (Mi , fi) planes, for the two possible case Mi < M2 (in orange) corresponding to 
SPSla, and also for the alternative solution with M2 < M\ (in red). is determined with an accuracy of about 
A/i ~ 15 GeV, but as already mentioned the sign of /1 remains undetermined. Moreover, in a most general MSSM, 
without any prior knowledge on gaugino and Higgsino mass relative values, the sole knowledge of three neutralino 
masses does not determine the relative hierarchy among Mi, M2, and So strictly speaking there is a six-fold 
ambiguity in this case, considering all possible ordering of these three parameters (see Appendix A). 
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FIG. 7: Constraints on Mi, /i (at Qewsb scale) in general MSSM for three neutralino mass measurements in gluino + squark 
cascade decay. 



We now consider the gaugino mass universality case with three neutralino mass input. The ambiguity on the relative 
magnitude of Mi, M2 and fi obtained in a general MSSM is of course resolved in this case since the hierarchy at the 
low scale is entirely determined from universal initial values of Mj. The resulting constraints on fx, tan (3 are illustrated 
in Fig. [8] where one observes that the /i < solution has disappeared (and the equivalent of Fig. [4] would show that 
only a smaller part of the red contour fj, > is surviving). However, only fx is much more constrained, while apart from 
a slightly more interesting lower bound, tan (3 remains essentially unconstrained for large tan (3. This is simply due 



to the sin 2/3 only dependence in all these relations, whatever the number of input neutralino masses. More precisely9 

tan/3>2.7, 350 GeV < fx < 372 GeV . (3.13) 

These results are consistent with general expectations, namely that the gaugino sector alone can hardly constrain tan j3 
at the LHC, even in the mSUGRA (gaugino universality) case, but these features are perhaps very simply illustrated 
here analytically. 




FIG. 8: Constraints on /x, tan/3 from gaugino universality and gluino plus three neutralino mass measurements in gluino + squark 
cascade decay. 



Concerning the inclusion of radiative corrections into this essentially tree-level determination, we have already 
mentioned the use of the DR parameters Syy, rhz in all these relations, which already incorporates a part of the cor- 
rections. The bulk of these DR corrections with respect to tree-level values come from standard model contributions, 
and supersymmetric contributions, though not negligible, are strictly speaking inducing some theoretical uncertainties 
if loop contributions from unknown sector are taken into account. Next, the corrections from the running to the pole 
neutralino masses remain moderate in the MSSM, so that one may neglect them in a first approximation. Actually, 
most of these radiative corrections can be incorporated in the form of corrections to the "bare" parameters n, Mi, 
M2 to very good approximation 13 811 . while additional corrections in the neutralino mass sector are smaller[38]. Our 
procedure to determine e.g. Mi, M2 or /i above remains thus correct, up to an implicit shift on these parameters, 
provided that one has some knowledge on these radiative corrections. In fact gluino mass corrections are dominantly 
due to squarks/quark, and quite similarly for neutralinos, so they can be evaluated within a reasonable approximation, 
without assuming more knowledge than the available cascade decay input 10 . Though it is not completely straightfor- 
ward to incorporate consistently all these corrections within these parameters, it is very similar to the kind of iterative 



NB the radiative corrections to Mi , M2 and fi could be used in principle to try to reduce the Mi < M2 versus M 2 < M 1 etc reconstruction 
ambiguities in the unconstrained MSSM. However these corrections are far too small to disentangle these ambiguities, at least for the prospected 
LHC neutralino mass accuracies. 



procedure required for a standard top-down calculation of MSSM spectrum |31|], as discussed in sub-section 2.3. 
is taken into account in the above illustrations. 



TABLE IV: Constraints on gaugino-Higgsino sector parameters (at Qewsb scale) from a standard top-down \ 2 MINUIT fit of 
neutralino and gluino masses under different universal or non-universal gaugino mass assumptions. 



Data & fitted parameter 


MIGRAD (68%C.L.) 


MINOS (68%C.L.) 


nominal 


(+ model assumptions) 






SPS la value 


m iV 1 > m JV 2 > m 9 


(convergent) 


(problems) 




(Non-universal M;; 








1-loopRGE; noff,R.C.) 








Mi 


-97-126 (from fit) 




108.8 


M 2 


— 181—381 (from fit) 




208.9 




200-1500 (scanned) 




340.6 


tan /3(Qewsb) 


1-50 (scanned) 




9.74 




(convergent) 


(convergent) 




Mi 


108.8 ±7.2 


108.8 ±7.2 


108.8 


M 2 


208.9 ±8.2 


208.9 ±8.2 


208.9 




340.6 ±11.6 


340.6 ±11.6 


340.6 


tan/3(Qevi/ss) 


9.74 (fixed) 


9.74 (fixed) 


9.74 


m ]V 1 ' m JV 2 . m 9 


(convergent) 


(problems with tan (3) 




(Universal M^Qgut); 








1-loopRGE; noiViR.C.) 








M 3 


603.8 ±13.5 


603.8 ±13.5 


603.8 




341.6 ±293 


341.6 ±293 


340.6 


tan P{Qewsb) 


9.73 ±52.6 


(not calculated) 


9.74 


m Ni > m N 2 ' m N4 ' W 9 


(convergent) 


(problems with tan /3) 




M 3 


603.8 ±13.4 


603.8 ±13.4 


603.8 




340.6 ±13 


340.6 ±13 


340.6 


tan(3(Q E wSB) 


10.0 ±15.3 


(not calculated) 


9.74 



In order to cross-check the inverse analytical determination above, we perform alternatively standard top-down \ 2 
minimization fits, using MINUIT. The best fit results, with or without gaugino mass universality assumptions and for 
different assumptions on /i and tan/3, are shown in TablellVlboth for two and three neutralino mass input. Actually, 
in the non-universal gaugino mass case, two neutralino masses are clearly not enough input to constrain the four 
parameters Mi, M2, [i and tan (3. Thus to compare as much as possible with the previous reconstruction results, in 
this case we also perform scans over /i and tan (3 (taking fi > 200 GeV to avoid the pole region discussed above). The 
corresponding results obtained for the "envelope" of best fit values, with confidence level domains for M\ and M2, 
are shown in the second column of Table |IV| (in the first two entries). The rather bad constraints obtained in this case 
are roughly consistent with the above analysis of the /i dependence of M2 and Mi illustrated in Fig. [2] and already 
discussed above. (Note that our simple top-down fit is not able to exhibit the Mi , M2 reconstruction ambiguities in 
the non-universal gaugino mass case, discussed in sub-sec. 3.1, 3.4 and Appendix A, though these ambiguities may be 
implicitly responsible for the large errors found for Mi, M2, when fi is varied in a wide range). Similar ambiguities 
are however more clearly exhibited when performing more sophisticated minimization, where typically extra "best 
fit solutions" appear lll7lP . In contrast, the determination of Mi, M2 is much improved if ji is better constrained, 



and accordingly there is a substantial improvement on the /i determination if the third neutralino mass input 
available, as shown in the Table. Those results from Table II VI are thus roughly consistent with the analytic behaviour 
illustrated in Figs. [3j [7] and [8] The symmetrical (MIGRAD) minimization, however, does not reflect very well the 
true sensitivity to some of the parameters, most notably for tan j3. Indeed the unsymmetrical non-linear MINOS 
minimization does not give any useful constraint for tan j3, even for universal gaugino mass assumptions. Note that in 
the latter case , are essentially Bino-like and Wino-like respectively, while m ^ is essentially Higgsino-like, 
which also explain the drastic improvement on \i accuracy from a third neutralino mass input. 

E. Reconstructing gaugino masses at GUT scale 
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FIG. 9: Bottom-up RGE of gaugino masses from EWSB to GUT scale, including error propagations (/i and tan ft are fixed to their 
SPS la values). 

Finally we perform a bottom-up RGE of the gaugino mass parameters in the non-universal case, in order to check 
eventually for their unification at a high GUT scale, following ref. lEoll . (NB the bottom-up RGE procedure for 
gaugino masses and other basic MSSM parameters from given values at EWSB scale is explained in Appendix B). 
This is first shown in Fig. [9] taking the best M^Qewsb) determination above (i.e. corresponding to three neutralino 
mass input), and in Fig. [10] taking the worst Mi(Q ewsb) determinations (i.e. a scenario with two neutralino mass 
input and no known constraints on /z). One can see that a very good check of GUT scale universality is possible as long 
as the initial M{ (Q ewsb ) accuracies are reasonable: in other words, the error dispersion from the gaugino mass RGE 
remains small (which is clearly explained from the fact that their RGE content only depend on gauge couplings at the 
one-loop level). As expected this is qualitatively consistent with the former results of ref. lEoll (though it appears that 
the Mi, Mi LHC accuracies considered at that time were slightly more optimistic than the ones we obtain here from 
our analysis). Note that the dispersion due to the bottom-up RGE can be much more important for other parameter 
sectors, in particular for the m,H u parameter in the Higgs sector (see the discussion in Appendix B). Nevertheless, the 
possibility of checking universality at GUT scale may become elusive even for gaugino masses in the extreme case 
where almost nothing is known on [i, such that the Mi, M<x low scale values errors are large. In this case only the 
Ms error remains under control, as illustrated in Fig.[l0] We also remark once more that in addition to performing a 




bottom-up RGE, mSUGRA gaugino universality could be checked efficiently from plots as illustrated in Fig. [2] (where 1 
it can also be eventually distinguished from other high scale SUS Y-breaking models like AMSB). 
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FIG. 10: Bottom-up RGE of gaugino masses from EWSB to GUT scale, including error propagations for the worst case (/i and 
tan (3 undetermined). 



IV. SQUARK/SLEPTON PARAMETER DETERMINATION (FIRST TWO GENERATIONS) 

We now consider a bottom-up reconstruction approach for the (first and second generation) squark and/or slepton 
sector parameter of the MSSM. Because of the negligible mixing for the first and second generations, as well as 
the simplified RGE structure of this sector, we can elaborate a rather general strategy to reconstruct the relevant 
parameters up to an (eventual) GUT scale. Similarly to the gaugino/Higgsino sector, we shall consider two different 
model assumptions: a general MSSM sfermion sector, or a more constrained scenario assuming universality of slepton 
and squark masses. 

A. General MSSM and a simple squark/slepton mass sum rule 

The masses of the sfermions b, t) and / which participate to the cascade decay in Eq. ( 12. It obey (at the EWSB 
scale) well-known relations (at tree-level) e.g for the up squark and selectron: 

1 2 

m l ± = m l L + (3 - ^s^ v )m 2 z cos2f3 

m i 2 = m l R ~ s w m % cos2 /3 

which are valid in the general MSSM . Similar relations hold for the other squark and slepton flavors. Eqs. (14. j} 
simply relate the physical masses to the soft breaking scalar masses via the additional D-terms. These relations are 



s^, in Eq. j4.lt and further equations below should be understood as the D_R-scheme parameter 
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particularly simple for the first two squark and slepton generations due to negligible mixing. The main idea is -?6" 
consider specific linear combination ("sum rules") in order to eliminate the tan/3 dependence: 

^"4, + i\ - \ s w) m l 2 = s w m l L + (5 ~ \ s w) m l R ■ (4-2) 
Then taking into account the available accuracy on the physical masses to Ui and mg 2 provides constraints on the 
MSSM soft-breaking scalar parameters ma L , m^ R independently of tan /3 values. So, even if tan /? is largely unde- 
termined (as is the case from the neutralino sector alone at LHC, illustrated in previous section), Eq. ( 14.21 ) gives a 
(model-independent) determination of the linear combination of basic parameters which will be roughly of the order 
of magnitude of the physical mass accuracies, i.e. a few percent for typical LHC prospects. More precisely, a straight- 
forward calculation from the experimental accuracy in Table Hill gives a ~ 2.3 % (~ 2 %) relative accuracy for the 
linear combination in Eq. (14.21 1. if we combine the and mg 2 errors linearly (quadratically), respectively. (NB we 
have neglected at the moment for simplicity the error on s^: the effect of Sy, uncertainties will be studied below). 
Another advantage of this linear combination is that, even in a general non-universal MSSM case, the RG evolution of 
the relevant parameters (ma L (Q), mg R (Q)) depends only on the gauge couplings gi and the gaugino masses Mj 12 . 
The linear combination Eq. ( 14.21 ) can thus be RG-evolved within a restricted set of input parameters solely determined 
from the gluino cascade, in order to obtain the soft scalar term values at the GUT scale, where one may check for 
eventual universality relations. Before doing this RG evolution, it is necessary to subtract out radiative corrections 
linking the running to the pole masses, which are not negligible for squarks, in the way discussed in section 2.3. More 
precisely we have 

M C e = "*ai(Q) + Ama, (Q, m q , m s , ...) (4.3) 

where these corrections are largely dominated by squark/gluino contributions at one-loop. For a typical mSUGRA 
scenario like SPSla we have Am Ul (Qewsb, mq, m§) ~ 19 GeV, which can be consistently subtracted out to define 
the running mass m Ul . Concerning radiative corrections linking the running to the pole slepton mass, they are gener- 
ally much smaller and we shall neglect them in our analysis. 

Next, it is a straightforward exercise to work out the RG evolution of the linear combination ( 14.21 ): 

d 



dt 



i 2 w m l L + (\ ~ \ s w) m l R 



2 dm k ,A 2 2 dm\ R 



where t = In Q and the standard RG evolution of the relevant parameters m\ , m 2 is used (which as mentioned only 
depend on gi, Mi). We also have: 

4t=^+^{l^4->*4) 

withs 2 v (<) = I g\ (t)/(§<7i (t) + g\ (<)) and the factor 3/5 is due to the standard normalization of the U(l)y coupling 
gi in the MSSM RGE. Eqs. (1441 ) and d43T ) take into account the RG evolution of s%j. (Note that the latter is not at all 
negligible since it is related to the running of gauge couplings, which change substantially from mz input values to 
GUT scale values). 

As already mentioned one important feature of this bottom-up RG evolution is that some of the low scale parameter 
uncertainties are amplified once evolved to a large scale, depending on the structure of RG beta functions for some of 
the relevant parameters: this is the case to some extent with the evolution of Eq. ( 14.21 as we will see later. 



This is only true at the one-loop level RGE, since at two-loop level practically all MSSM parameters enter the RG evolution of squark and slepton 
mass terms [26l |27jl . We shall discuss below how the inclusion of two-loop RGE effects affect our results, but we anticipate that neglecting these 
higher loop effects do not change drastically the obtained constraints. 



B. Constrained MSSM with squark, slepton universality z:> 

If we now assume squark and slepton mass universality at the GUT scale, Eq. (14. A\ immediately determines 

ml' (Qgut): 

s 2 W ™l L +{\~ \*wWe R (Qgut) = | K'<) 2 (4.6) 

where the gauge coupling universality relation at the GUT scale: Sjy(t = IuQgut) = 3/8 has been used, mj 1 
indicates that we only assume universality for the (first two generation) squark and slepton sector at this stage, i.e. not 
necessarily for the third generation sfermions, nor for Higgs scalar terms like in mSUGRA models. 



C. Explicit reconstruction test for the SPSla input 

We now determine explicit constraints on the squark and/or slepton sector parameters from the specific SPS 1 a blind 
input with the expected accuracy on the masses of and Ir from Table HIT1 As already mentioned in section 2, there 
is at present no way to tag the charge and flavor of the relevant (first two generations) squark at LHC. Accordingly 
the resulting mass accuracies of say, an up or down squark, are assumed to be identical [4], so that there is no need 
to combine their errors in a statistically elaborated manner, and we thus assume that it is sufficient for our purpose to 
take the average of two (identical) errors in our analysis. By combining thus the accuracies on the measured u\ and 
(=2 masses in Table HTT1 we obtain for SPSla from Eqs ( 14. 2K (14. 6t : 

84 (86) GeV < m q l < 116 (112) GeV . (4.7) 

Note that the first limits are for linearly combined mass uncertainties (while those in parenthesis are for quadratically 
combined mass uncertainties). We emphasize that the bounds in Eq. (14.71 ) are independent of tan (3 values. However, 
there is a rather important amplification of the initial low scale uncertainty ~ 2% due to error propagation in the 
bottom-up RG evolution, and because of the additional terms proportional to ds^/dt in Eq. d4.41 i. To better trace 
the origin of the resulting mo uncertainties, it is illustrative to consider independently the squark and slepton mass 
uncertainties: this gives 

86.6 GeV <m^(Am fil ) < 112 GeV (4.8) 



and 



97.1 GeV < mg''(Amg a ) < 103 GeV (4.9) 



respectively. Thus the final uncertainty on mo is largely dominated by the initial accuracy. Actually, the latter 
bounds on ra^ 1 were calculated while fixing the gaugino mass terms Mj. But since the RG evolution of mn L , mg R 
depends on gaugino masses, the mj' 1 accuracy should be sensitive to Mj uncertainties, mainly those of M3 which 



are enhanced in the m„ L RGE by the strong coupling: dm 2 UL /dt ~ — ^agM 2 |[26t 12711 . This leads to an important 



amplification of m^' final uncertainty due to the AM3 ~ 7 GeV uncertainty (although the latter effect is damped 
somehow by sfy ~ .238 in the first term of the RHS of Eq. ( 14.41 ). while other terms in Eq. ( 14.41 ) are not much sensitive 
to Mj uncertainties). One thus obtains, in the conservative case of uncorrected and linearly combined errors, a 
maximal uncertainty on m^' of about ± ~ 32 GeV (respectively ± ~ 22 GeV when combining errors quadratically) 
instead of the bounds shown in Eqs. (14.7b — (14.9b 13 . However, the linear combination Eq. ( 14.21 ) does not use the 



Note also that the dominant uncertainties from mj L1 and rrig (i.e. M3) have the opposite (anti-correlated) effect: from the RG structure, increasing 
(resp. decreasing) makes mj' 1 to increase (resp. decrease), while the opposite behaviour is obtained from M3 ± AM3. 



full information from the two independent mass relations in Eq. (14.11 ): we shall illustrate below how this additional 
information improves rather substantially the limits on TOq''. 

Another potentially interesting question is whether one can derive at the same time any useful limits on tan j3, once 
using the complete information from both squark and slepton masses. At first sight one may naively expect to obtain 
some upper bounds on tan j3, since the relations d4.lt are sensitive to cos 2/3. However a simple estimate immediately 
indicates that interesting tan (3 upper bounds from this squark, slepton sector are hardly expected for the given LHC 
e.2,U\ mass accuracies: in fact for the SPS la point with tan (3(Qewsb) ~ 9.74, cos 2/3 ~ —0.979 i.e. very close to 
— 1, so that one would need at least an accuracy < 0.02 on | cos 2/3| to put useful upper limits on tan /3. In contrast, a 
simple calculation of uncertainties from both Eqs. ( 14. U gives: 

A|cos2/3| - 0.6 (4.10) 

even in the optimistic case where we combine the mg I , mg 2 uncertainties quadratically, and neglect the errors on sfy. 
The above estimate, however, does not take into account other possible constraints on tan (3, that may come from other 
sectors, or from theoretical consistency. For instance the obvious constraint: | cos 2(3 \ < 1 puts additional limits on 
2 — *"| via Eqs. ( 14. It , so indirectly on m^' . Furthermore, for uiq' 1 ~ 100 GeV it is easily checked that tan/3 



m e R - m e 2 

cannot be larger than tan (3 < 35 — 36, since beyond this value the lightest stau mass becomes tachyonic due to 
the large stau mixing fh T ./i tan (3 term. Moreover, the lightest Higgs mass becomes inconsistent (or too low) for small 
tan (3 < 2.2 approximately, and for low tan/3 values the LEP lower bounds on m^|45|] actually put a tighter constraint 
tan (3 > 8 — 9. Yet the latter limits are theoretical and model-dependent (or experimental and model-dependent in the 



case of mh bounds), and specific of the SPSla benchmark ll34ll . which was chosen on purpose to satisfy the present 



experimental constraints. If we push m 1 ^ 1 sufficiently above the SPSla central values, the upper bound tan/3 < 35 
from tachyonic fi is easily evaded, though in this case one should also take into account the lower bound > 104 



GeV from LEP limits 14611 . (For example, for tuq ~ 200 GeV tan/3 ~ 50 is not excluded by tachyonic stau, while 



the gluino, squarks cascade decays would not be drastically different from the SPS 1 a one). Similarly, the tan (3 lower 



bound due to LEP m/i lower limits could easily be evaded in an unconstrained MSSM[4?]. We will thus not apply 
such direct (or indirect) experimental limits which are much dependent on the specific SPSla choice, since our main 
aim is to present a reconstruction strategy expected to be valid beyond this particular benchmark choice. 

Accordingly a question that we examine in some detail next is whether the sole , mg 2 mass measurements could 
put some extra model-independent limits on tan/3. From the previous estimate it appears that to obtain stringent such 
experimental constraints on tan j3, one would require an accuracy about an order of magnitude better on the squark 
and slepton masses than the one prospected at the LHC. Incidentally this is roughly the accuracy expected at the ILC 
(though only for the sleptons), where both Ir and li masses could be measured at the per mille leveljsi Hi - However, 
a detailed ILC analysis is beyond the scope of the present paper and left for future work. 

We anticipate that (model-independent) limits on tan (3 will be indeed absent, or very marginal if using solely the 
first two generation squark and slepton mass accuracies. This is consistent with general prospects flS]- We found 
however useful to examine this issue in some detail, since our construction is not limited to the LHC mass accuracies 
here considered: thus tracing analytically the sensitivity on parameters can help to understand better what determines 
the constraints in a more elaborated analysis. 

To begin, there is a subtlety that is not taken into account in the above crude estimate of error combination in 
Eq. d4.10| ), such that it may underestimate the tan (3 sensitivity: as emphasized previously the RG equations for 
ht-iil , m e R only depend on gauge couplings and gaugino masses, so in particular do not depend on tan (3 (at one-loop 
level). But there is in fact an indirect dependence on tan (3 even at one-loop (though very moderate): it originates from 
the boundary conditions on the RGE, namely the initial values of gauge couplings gi in the DR scheme, as well as their 
values at the GUT scale (if gauge unification is imposed), depend slightly on tan (3 through radiative corrections (see 



e.g. ref. 13811 ). Although these effects are strictly speaking small higher order corrections on gi(mz), giimaur) val- 



ues, they are enhanced when running parameters from low to high GUT scale over more than 13 orders of magnitude. 



Indeed most realistic calculations of MSSM spectra do take into account this dependence consistently! 28, 29, 30, 3 iff 
Also, these values are notoriously different if using the RGE at the full two-loop level, or in a one-loop approxima- 
tion, as will be illustrated below. (This is somewhat similar to the impact of precise initial gauge couplings and RGE 
approximations on the GUT scale value ~ 2 x 10 16 GeV). Moreover, the low energy EWSB scale at which all soft 
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FIG. 11: ma L (QEWSB) as function of tan /3 for SPS la in one-loop RGE (top plots) and two-loop RGE (bottom plots). Red dashed 
line: one-loop RGE, default Qewsb = (mj l mf 2 ) 1//2 . Black line: fixed Qewsb = 468 GeV. Orange dotted line: two-loop RGE, 
defeat Qewsb = (m^m^) 1 ^ 2 . Blue line: fixed Qewsb = 468 GeV. 



parameters are evaluated also depend in principle on tan (3 values, and more generally it is not fixed by first principles: 
in most analysis the default EWSB scale is often fixed to Qewsb ~ i 771 ^ 171 ^) 1 ^ 2 (which is know to minimize the 
scale dependence 148J of the MSSM one-loop effective scalar potential). For the SPSla point this is a well-defined 
value, but when relying only on the gluino cascade decay data, the stop masses are not assumed to be know in a general 
MSSM case. Thus we can in principle vary this EWSB scale, which can affect the final value of soft scalar masses, 
since this scale determines the endpoint of the RG evolution. In our calculation we use either a fixed value (close to 
the true SPSla for definiteness), or adopt the above default value in universality cases. In fact m^ L is expected to 
be specially sensitive to such variations, due to the large coefficient ~ as-Mf within its RGE, which makes it to run 
much faster than m e - B . (E.g. for SPSla m^ L (Qgut) = 100 GeV evolves to rriu L {Qewsb) ~ 560 GeV.) 

For illustration in Fig. QT| we vary rriu L (Q ews b ) as function of tan (3, at one- and two-loop RGE order, both with 
fixed and default EWSB scale, using SuSpect 2.41. The variation with tan/3 is less than 2 GeV for m UL (and we 
checked that it is completely negligible for m eR ) for the whole tan (3 theoretically allowed range. This is accordingly 
below the prospected LHC accuracy on m^, and thus rather negligible for our analysis. In contrast, there is a large 
difference between one and two-loop RGE, but this is usually the case for the whole MSSM spectrum, as it is well 
known. This illustrates that for such reconstruction (or in fact any other reconstruction methods) one should be 
careful to be consistent with the RGE approximation used. Concerning now the small tan (3 dependence it should 
be kept in mind that any more elaborated fit (as could be performed from a x 2 minimization using MINUIT) will be 
eventually sensitive to such effects, since these cannot be easily "switched off" from the fitting procedure. Moreover, 
such indirect dependence on tan (3 will be relevant anyway once reaching a better accuracy on slepton masses as is 
prospected at the ILC. 



All these features are examined more quantitatively from a systematic scan over parameters, where we took into 



account the propagation of rrif ll , mg 2 uncertainties in both relations d4.ll ) after evolving from low to high scale, afi§ 
back. More precisely we scanned over allowed mass bounds, considering the two mass uncertainties as independent 
(and uncorrelated) for simplicity. We also took into account any additional small dependence e.g. on tan j3 such as the 
one illustrated in Fig.QT| The (mg' J ,tan/3) distributions obtained are shown first in Fig.Q~2](the scan was performed 
with 3000 uniformly distributed random numbers). Comparing with Eq. j4.7b one notes that the lowest values of 




FIG. 12: Constraints on mo = mj' 1 (Qgut) and tan/3 from first two generation squark and slepton in gluino cascade decays 
after bottom-up RG evolution. The orange (respectively magenta) points correspond to constraints obtained from the squark ma 1 
(respectively mg 2 ) relations in Eqs. l l4.lt . 



84 < too < 94 GeV have been excluded simply from the constraint | cos2/3| < 1. However, arbitrary large values 
of tan/3 are possible (provided that a sufficiently large number of scan points are taken), confirming the simple error 
estimate above in Eq. ( |4.10t . The fact that large values of tan (3 are very few in Fig. [12] while low tan (3 values appear 
very much favored, is actually an artifact of the (uniform) scanning procedure where we basically scanned over the 
mg I , TOg 2 masses, which determine cos 2/3 via Eq. d4.U rather than tan/3. So the distribution of points in Fig. Q~2] 
is simply resulting from the transformation from cos 2/3 to tan/3 and has not much statistical meaning as "most 
likely" values of tan/3. As already mentioned we shall refrain to enter into a fully realistic statistical treatment of 
uncertainties, which would require to take into account lfl7ll . among other things, the non-trivial correlations implied 
by the cascade decay mass measurements. It is yet tempting in the present case to proceed one step further with a 
little more elaborated analysis. We thus perform a different scan, where instead of uniformly distributed "flat prior" 
random numbers, we exploit simple trigonometric relations to match more faithfully the true tan (3 distribution. In 
addition we start with Gaussian-distributed random numbers (assuming thus that the two independent errors are 
purely statistical, which is actually not really correctly 0]). We then calculate a "theoretical" x 2 (where Xmin = 
trivially since we use the correct central values of the masses) just to obtain well-defined confidence levels (C.L.) for 
the (joint) estimation of the two parameters m^' 1 and tan/3, assumed to be independent and uncorrelated. The result 
of this Gaussian scan is shown in Fig. Qj] The domains obtained from the TOg 2 and to^ relations in Eq. ( 14. U . at 
one-a level (more precisely 68% C.L.) are shown in black, and the additional points allowed at two-er (95%) C.L. 
are shown in indigo and orange respectively. The distribution of those allowed points is now essentially uniform in 
tan (3 (though now the apparent concentration for relatively large tan /3 is purely an artifact of the logarithmic scale), 
in contrast with the previous plot in Fig. Q~2] More importantly, tan (3 is constrained by lower and upper bounds, but 
these essentially originate from the theoretical (model-dependent) constraints, indicated in red for large tan (3 which 
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correspond to a tachyon f\. (NB the white zone for tan (3 < 2.2 corresponds to a very light m/J. More precisely we 
observe that even at the 68% C.L., values of tan/3 up to the theoretical SPSla constraint tan /3 < 35 are not excluded. 
Moreover one sees that the 95% C.L. additionnal points (indigo and orange) have only the effect of slightly enlarging 
the mj' 1 determination but not much influence on the tan/3 determination. Thus a more realistic treatment of errors 
simply confirms our above crude estimate, indicating that no interesting model-independent upper limits on tan (3 can 
be derived from the prospected LHC squark and slepton masses accuracies. 

It is instructive to compare at this stage those results with a more standard top-down fit. We thus used MINUIT to 
perform a standard \ 2 minimization starting from a MSSM model with universal GUT mg' 1 and tan (3 free parameters. 
First we fix the gaugino masses Mi, M2, M3 to their SPSla values and then fit this model to the "data" consisting 
solely of 1 TO e 2 ■ This two-parameter fit result is shown in Table[V]for different MINUIT minimization options, at 



the one-er level. In fact the simpler MIGRAD ll 1H (symmetric error) minimization did converge, giving apparently a 
(marginal) upper tan (3 limit at the one-a level. But it is not a very useful tan (3 bound, being above the theoretical 
tan/3 limits in this SPSla case. Moreover, the more elaborated MINOS minimization, taking into account properly 



unsymmetrical errors and non-linearities 11 111 , did not find any tan/3 upper (nor lower) limits. These results are thus 
qualitatively consistent with our more naive analysis above. Note that the mo bounds are also roughly consistent with 
our previous results. Now there are several reasons not to trust even the marginal upper bounds found for tan (3: though 
it is difficult to trace the very details of the minimization steps, the fit is much probably indirectly influenced by the f 1 
becoming very small (and ultimately tachyonic) for tan (3 > 35 — 36. More precisely, above those values, SuSpect still 
gives m Ul , nig 2 output (unless explicit warning flags are switched on), but these are no longer very reliable (because 



the iterations needed to calculate a convergent spectrum are stopped in this case 113 111 ). In particular, there is an abrupt 
change of m UL once tan (3 > 36, see Fig. [TTJ (though the overall variation remains reasonable). Thus, comparing these 
minimization results and Fig. [13]with the simple estimate Eq. ( 14.1 OK we can infer that the marginal upper bound on 
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TABLE V: Constraints on mo, tan/? obtained from a standard top-down \ 2 nt with different model and input assumptions and 
different level of MINUIT minimizations. 



Data & fitted parameter 
(+ model assumptions) 


MIGRAD 
minimization (68%C.L.) 


MINOS 
minimization (68%C.L.) 


nominal 
SPS la value 


(1-loop RGE + nogR.C.) 

a I 

tan /3 


(convergent) 
99.98 ±8.8 
9.41 ± 29.8 


tan /3 exceed limits 

99.98 ± 8.7 
q ^ +no limits 
-no limits 


100 
9.74 


mg 2 ,m fil + 
(1-loop RGE + noqR.C) 

n,l 

m o 
m 1 / 2 
tan /3 


(convergent) 
99.96 ±7.9 
250 ±4.3 
10. ± 29.8 


tan j3 exceed limits 
qq q +io 

yyy -10.5 

250 ±4.3 

J0 +29.8 

-no limits 


100 

250 
9.74 



tan (3 is principally determined by the indirect small higher order tan /3 dependences (like typically the dependence 
ma L (tan/3) as illustrated in Fig. QT| originating mainly from the slightly varying EWSB scale, rather than directly 
from the mj,, mg 2 accuracies). We have crosschecked this by redoing similar fits with a constant EWSB scale, or 
with two-loop RGE, which both have the effect of smoothing somehow the variation of rriu L (tan (3) near the transition 
to tachyonic f\ for tan (3 ~ 36, and the tan (3 upper bound in TablefVltends to increase (or even to disappear with non 
convergent minimizations). We have also further checked this by increasing progressively the mass accuracies: while 
the corresponding mo bounds decrease, following the expected statistical behaviour, for tan (3 one obtains either non 
convergent minimizations, or extra odd solutions far from the SPS la values. There are anyway no improvements on 
tan (3 bounds until a substantial decrease of these experimental errors is set (about an order of magnitude smaller than 
the LHC accuracies of Table [TTTb- 

For completeness we performed another minimization, taking into account in addition the neutralino and gluino 
masses with a three-parameter fit of mo, TO1/2, tan/3. The situation does not improve much as concerns tan (3 limits, 
as illustrated by the corresponding results in TablelVl In this case the (unsymmetrical) MINOS mo bounds are slightly 
worse than for the two-parameter fit, which is expected since now A/3 is not fixed. Accordingly the M3 accuracy 
propagates to the mo determination via the RGE, as explained above. On the other hand the accuracy on is very 
good. Finally we also performed a similar fit for a benchmark point like SPS la except mo = 200 GeV, in which case 
the transition to tachyonic happens only for very large tan (3 > 50, having thus potentially less influence on the 
fit. Assuming the same mass accuracies than for the true SPS la point we obtain, for a two-parameter (respectively 
three-parameter) fit: tan/3 = 9.93 ± 52.2 (respectively tan j3 = 10.1 ± 50.7), at one-cr level, while errors for mo and 
m 1 / 2 are very comparable to the SPS la ones. So this confirm the above analysis and indicates that the tan (3 upper 
bounds are essentially inexistent. 

This also illustrates that a "global" top-down fit, whatever elaborated with MINUIT algorithms, may be fooled 
and lead to misleading conclusions due to extra parameter dependences which originate from theoretical approxi- 
mation artifacts (e.g. here the choice of EWSB scale, RGE approximation, etc). In the present case, one could of 
course easily avoid such problems by simply adding protections within the minimization procedure, but that would 
amount to put explicitly the (model-dependent) upper bound on tan (3 due to tachyonic (or more generally too light) f 1 . 

We finally examine two questions related to theoretical uncertainties that are relevant to the above analysis. First, 



as already mentioned in deriving the Mq' 1 constraints in Eqs. ( 14.7114. 8b we had neglected for simplicity the errors o'A 
Syy\ the latter are actually not quite negligible, since even a small uncertainty can in principle affect our determination 
from Eq. ( 14.21 ). However, the bulk of radiative correction contributions to the s^y parameter, in the DR scheme, 
originates from standard model and are thus predictable in our framework. Additional supersymmetric contributions 
are not negligible either ID8I1 . but we checked that varying all MSSM parameters and sparticle masses form the SPSla 
to models with arbitrary MSSM values produces a variation of of about 0.6% only. One may probably push 
parameters to extreme values to find a slightly larger variation, so we conservatively consider a 1% uncertainty on 
Sw. The impact on mo determination is an (upper) shift by about ~ 3 GeV with respect to the numbers quoted in 
Eqs. ( 14. 7114.8b . Note that the correct dependence of mo and tan (3 upon is automatically taken into account in the 
MINUIT fit results in Table |V] 

Another potential question is that the determination of m^ 1 via Eq. (14.21 i together with RGE, depends only on the 
cascade masses if restricting the RGE to one-loop order. At the two-loop order, practically all other MSSM parameters 
are entering the RGE. Nevertheless, it is possible to study the impact of this uncertainty by assuming, within the two- 
loop RGE level, simple (e.g. universal) relations for the unknown parameters, and to redo our analysis: while central 
values are evidently shifted, the impact on error propagation is rather negligible, with minor quantitative changes 
on e.g. the obtained m^' 1 constraints. This is also consistent with some comparisons we made of two-loop versus 
one-loop RGE fit results using MINUIT. 

We thus conclude that tan (3 is essentially unconstrained from the data we used at this stage, which is not much a 
surprise and consistent with general expectations on LHC prospectsji, 22]. Model-independent constraints on tan j3, 
though moderate, may be obtained however from other sectors as we shall see in next sections, either considering the 
information from bottom squarks in the cascade, or adding the lightest Higgs mass measurements. 



V. THIRD GENERATION SQUARK PARAMETER DETERMINATION 

We will now consider the possible determination of some of the third generation squark parameters. As discussed 
in section 2, the sbottoms enter the gluino cascade and both mass eigenstates may be measured to some extent (though 
with less accuracy for the heaviest &2)0 Etl- We will examine what additional information they provide, both in 
unconstrained MSSM or assuming universality relations, as in previous analyses. 

A. Scenario S5: constraints from the sbottom masses in non-universal MSSM 

So far, we have obtained from the gluino cascade useful constraints on the gaugino/Higgsino parameters and m^ 1 , 
while tan (3 is very poorly determined. The sbottom sector may provide information on the missing third generation 
soft scalar terms, or eventually on tan (3. To set up signs and other conventions, we recall the sbottom mass squared 
matrix: 



M 2 _ | m Q3L + m l + i~\ + |s^)m|cos2/3 m b (A b - /j, tan (3) 



2 3L T " 

rrib(Ab — ji tan (3) m b R + m b ~ § m z cos 2/3 



where again all parameters are implicitly understood to be in the DR scheme. It is immediate that the trilinear term Ab 
will be very badly determined, being largely suppressed by the bottom mass mb, with ^(Qewsb) ~ 3 GeV e.g. in 
SPSla, and similar values in more general SUSY scenarios. However, this also means that we can practically neglect 
Ab- In this case the sbottom masses determine, to a good approximation, the third generation parameters toq3 l , nib R 
in a simple way directly from (the trace of) ( 15.11 ): 



2 

mi 



m l + m b 2 + 2 C ° S 2/3 ~ ' h ^ S = Q3l + bR ' 



Now, taking one of the mass eigenvalue (say m^) as a given function of the other parameters fj,, tan/3 etc, we cilfi 
solve a simple relation for mq 3i , m bR : 



m Q3 L (b R ) 
D 



S+(-)D 



1/2 



-Y 



( m L _ m i - 2m b X b )(m? b2 - m? + 2m b X b ) 



1/2 



(5.3) 



where Y = {— | + |s^)to^cos2/3 and = — /itan/3. There is in fact a twofold ambiguity, namely a 
second solution where mQ 3L <-> in Eq. ( 15.31 ), since it is a second order equation. For SPSla we know that 
m Q3L (Qewsb) > m bR (QEWSB) resolves this ambiguity, but for unconstrained MSSM one has to deal in principle 
with the two possibilities in the absence of further information other than the sbottom masses. In addition one also 
expects radiative corrections to this tree-level determination which might spoil its simplicity. Radiative corrections 
to sbottom masses are not negligible in general, but to a very good approximation they are largely dominated by 
gluino and squark contributions, as already discussed in sub-section 2.3. So we have a reasonable knowledge of these 
corrections, such that the running masses in the DR scheme involved in Eqs ( 15.3b may be obtained from the pole 
masses after appropriate subtractions. 



B. Explicit reconstruction: SPSla test case 

The resulting rriQ3 L , m bn determination for the SPSla test case with corresponding input mass accuracies is illus- 
trated in Fig. [14] We plot the allowed range for rriQ3 L , m bft obtained from a (uniform) scan over mass accuracies 
from Table Hill and the other relevant parameters, as indicated. Different cases are illustrated when varying tan/? and 
H within the previously obtained limits, i.e. // ~ 360 ± 200, or more constrained situations, A/i ± 10 (the latter 
corresponding roughly to the case where three neutralino masses can be measured), and 3 < tan (3 < 35 (using the 
theoretical upper bound for SPSla). Since A b is essentially unknown at this stage, we vary it widely in the range: —1 
TeV to 1 TeV. As one can see, the determination of rriQ3 L , m ba is quite accurate even in a non-universal MSSM case, 
and for largely unknown A b . This is clearly explained by the strong suppression of A b , which has thus practically no 
influence on the scalar mass parameter determination. Moreover, the variation of \x, and even the large variation of 
tan /3, have no strong influence either, which is explained by the fact that they are both very suppressed also by m b 
as compared to other parameters. To examine the impact of radiative corrections, we did a similar scan but neglecting 
on purpose the radiative corrections, i.e. using as central values the pole sbottom masses. As a result the rriQ3 L , m bR 
central values are substantially shifted with respect to nominal SPSla values, but the accuracy obtained on mQ3 L , 
m bft is very similar. 

Next, in Fig. [15] we compare the previous results obtained from a simple scan with uniformly distributed random 
number with those from a scan made with Gaussian-distributed random numbers, defining a \ 2 domain similarly to 
the analysis in section 4. One can see that the corresponding one-cr limits show the expected statistical behaviour, 
having a more regularly-shaped contour than the one in the uniform case, but with no drastic changes in the limits 
obtained on niQ3 L , m bR . 

C. Scenario S' 5 : constraints from sbottoms in MSSM with sfermion mass universality 

Prospects for the third generation scalar terms are better in a more constrained MSSM i.e. with universality relations 
assumed in the scalar sector. In this case, one has: 



ml' = niQ3 L (Qgut) = m bR (Q GUT ) 



(5.4) 
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FIG. 14: Constraints in general MSSM on rriQ SL , m bR (at Qewsb scale) from m^m^ measurements in the gluino cascade 
decays. Radiative corrections linking the pole to running masses have been subtracted out. Experimental errors on sbottom masses 
are taken from Table Hill Blue region: fixed nominal values of tan/3 ~ 9.74, fj, ~ 357 GeV (and At = 0). Red region: 
3 < tan (3 < 35, A^ ~ 10 GeV, -lOOGeV < A b < WOGeV; green region: 3 < tan/3 < 50, Afi ~ 200 GeV, -ITeV < A b < 
ITeV. Nominal SPSla values are rriQ3 L = 503 GeV, m bR = 525 GeV (at one-loop RGE level). 




FIG. 15: Comparison of uniform (in green) versus Gaussian (in red) scan with one-cr (68% C.L.) level contours. The variation of 
extra parameters correspond to the less determined range in Fig. [14] i.e.: 3 < tan j3 < 50, A/i ~ 200 GeV, —ITeV < A b < ITeV. 



Then, using the above results from section 4 on jtIq limits, one can determine from a top-down RGE the corresponding 
accuracy obtained on m Q3L (Qewsb) and m bR (Qewsb)- 

m Q3L (Q E wsB)~ 498 ± 1.2 ±7GeV, m bR (Q EWSB ) ~ 521 ± 1.8 ± 6GeV , (5.5) 

where the first errors correspond to a variation of mo for fixed M3, while the second additional errors take into 
account the m~ g ~ M3 uncertainty (which dominates the final uncertainties). Alternatively, one may use here the 
limits on toq obtained from the \ 2 minimization in Table [V] which would gives results roughly comparable to the 



bounds in Eq. ( 15.5b (except that it is somewhat more difficult to disentangle the effect of the M3 uncertainties from trie" 
fit). Comparing with the unconstrained MSSM determination from the sbottom masses in Figs.[l4][l5]one can see a 
definite improvement on m,Q3 L (Qewsb) and mi )R (QEwSB) accuracies by about a factor two. This may be sufficient 
to resolve the twofold ambiguity between m,Q3 L and m& H discussed above in the non-universal sfermion case. 

Perhaps more interestingly, one may expect to derive an independent determination of tan (3 in this sfermion uni- 
versality scenario, by combining all information from too and the sbottom masses. More precisely, turning the other 
way round Eqs. (15.3b one can determine very simply Xb = Ab — [i tan (3 in terms of the know parameters, as follows: 



2 m b X b = - 



1/2 

(5.6) 



(Actually Eq. ( 15.6b comes from a second order equation, so in principle there is again a twofold ambiguity: Im^X^ = 
±[- ■ -] 1 ' 2 . But within our sign convention, if fi > 0, X b < necessarily, even for general MSSM, unless tan [3 < 
A b , which is not a very common situation except perhaps for very small tan (3 values.) 14 
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FIG. 16: Constraints on tan /3 with uniform scan from , accuracies. Ab is assumed to be almost undetermined (Ab = 
A(SPSla) ± 1 TeV). In blue : Am-^ = 7.9 GeV. In orange: Am 6 , = 16 GeV. 



The resulting determination of X b and tan j3 from Eq. ( 15.6b is first illustrated in Fig. Q~6j while Xb is directly 
determined from Eq. ( 15.6b , for tan (3 one has to take into account the additional uncertainties in the determination of /x 
(and Ab, which is assumed to be essentially unknown at this stage). However it is clear that this has a moderate impact 
on the tan/3 determination, since even a large variation —1 TeV < Ab < 1 TeV has a moderate effect, ~ AAh/fi 
(unless if would be very small). So even at this stage where Ab is completely undetermined, reasonable constraints: 
tan /3 < 27 — 28 are obtained, as shown by the blue contour. Moreover, even when increasing the uncertainty by 
a factor of two, i.e. Atoj ~ 16 GeV, one still obtains some reasonable constraints on tan/3 < 30, see Fig. [16] (But 
the upper bound on tan (3 disappear if the 62 is not measured at all). Like for the previous section analysis, we also 
show the difference between uniform and Gaussian scans, with corresponding one- and two-tr contours in Fig. [17] For 
completeness the joint (correlated) determination of (tan /3, Xb) is also shown in Fig. [18] 

We thus observe that, provided the two sbottom masses can be measured with this accuracy at the LHC, the bounds 
on tan (3 may start to be rather interesting, at least at the one-cr level. These results are also confirmed by a x 2 mini- 
mization fit using MINUIT, with more or less comparable bounds obtained, as shown in TablelVIl We finally mention 



Note also that tan (3 enters both sides of Eq. )5.6) via the D-term Y = (— ^ + \s^ v )vn? z cos 2/3. Thus we have to iterate on tan (3, and this is 
converging fastly. 



35 



liy[GcV] 




FIG. 17: Constraints on tan f3 from Gaussian scan: indigo: one-tr (68% C.L.); in orange: two-a (95% C.L.). 



X^GeV] 
n 



-3000 



-6000 



-9000 



-12000 



FIG. 18: Constraints on tan X b = At — /itan f3 from Gaussian scan: indigo: one-cr (68% C.L.); orange: two-a (95% C.L.). 



TABLE VI: Constraints on mo, tan (3 obtained from a standard top-down % 2 nt °f squark +sbottom masses and neutralino/gluino 
masses, for different options of MINUIT minimizations. 



Data & fitted parameter 
(+ model assumptions) 


MIGRAD 
minimization (68%C.L.) 


MINOS 
minimization (68%C.L.) 


nominal 
SPS la value 


mg 2 ,m fil , m- bl b2 + 
(1-loopRGE + nogR.C) 

q,l 

tan f3(mz) 


(convergent) 
99.96 ±9.9 
250 ±3.6 
9.98 ± 13.6 


tan f3 exceed limits 
qq q +io 

-10.4 

250 ±3.6 

9.98 +13 v . t 
-no limits 


100 
250 
10 




that squark and slepton universality relations as in Eq. ( 15.4b is yet not enough to have any interesting determination of 
At, which is quite obvious from its very suppressed contribution to the mass matrix d5.lt . One thus needs a different 



strategy to determine, if possible, the remaining trilinear couplings. The prospects for the latter would be certainly" 
much better if one could measure the stop masses in similar cascades. This is not the case for the SPSla, but we will 
examine in the next section how the Higgs sector parameters can help to determine the stop sector trilinear couplings. 



VI. BOTTOM-UP RECONSTRUCTION OF HIGGS SECTOR PARAMETERS 

In the Higgs sector, the most relevant parameters (at tree-level) are the scalar mass terms m# u , m H d and tan (3. 
These may be replaced equivalently, using the EWSB conditions, by /i, the pseudoscalar mass ttia, and tan (3. More 
precisely, the EWSB minimization equations dVuig g s/dHd = 0; dVm ggs /dH u = (where VHi ggs designates the 
MSSM scalar potential) can be solved for /j 2 and B/.r. 

fj, 2 = - [tan 2[3(m 2 Hu tan [3 - m 2 Hd cot (3) - to|] 

Bn = i sin 2/3 [m 2 Hu + m 2 Hd + 2fi 2 ] (6.1) 



which eventually includes radiative corrections to the scalar potential in the form of tadpole corrections t u , t d 1381,14911 : 

m 2 Hu — » m 2 Hu - t u /v u and m 2 Hd — * m 2 Hd - t d /v d ■ (6.2) 
Next, the running toa(Q) is defined as: 

fh 2 A {Q) = m 2 Hd (Q) + m 2 Hu {Q) + 2/z 2 (Q) (6.3) 
(where as previously all parameters are implicitly understood to be in the DR scheme). The A pole mass is then 



related to m^Q) via additional radiative corrections ll38ll40ll . It is well-known, however, that in the MSSM the scalar 
top sector contributes largely to radiative corrections to rah and niA- Thus even if some parameters are available from 
other sector analysis (though not precisely for tan (3 as we have seen), only an independent measurement of ttia (or 
alternatively an information on the stop sector, on X t = A t — [ij tan/3) could give a more useful information on 
the remaining parameters. For the determination of these Higgs parameter sector in a general (non-universal) MSSM 
case, the prospects are thus not very good if the available data are similar to those of the SPSla benchmark, even if the 
lightest Higgs mass could be determined with a good accuracy. In a constrained MSSM with complete high scale 
universality of scalar soft terms, prospects are possibly better. We can then use our previous analysis with in addition: 

ml' = m Hu (Qgut) = m Hd {Qgut) (6.4) 

to determine fi via Eqs. d6.1| ) as well as the running rriA value in Eq. (16.31 1. From Eq. d6.lt one has very roughly 
for tan/3 ^> 1: /j 2 ~ ~ m \ which is largely insensitive to tan/3. (NB in our analysis we use of course the 
complete expressions from Eqs. d6.1| l). Thus the rather good determination of too from the squark and slepton sector, 
as analyzed in previous section, together with universality assumptions d6.4t . give a more precise determination of 
[AQewsb)- However, tojj u has a strong sensitivity to the trilinear stop coupling A t through its RGE via terms 
ex y 2 A 2 , which restricts a very good determination of /i. Typically for 4 < tan j3 < 50 and a very moderate variation 
of A t , AA t ~ 100 GeV, we obtain: 

330GeV <ij(Qewsb) < 360GeV (6.5) 

which is comparable to the accuracy obtained in a general MSSM from neutralino mass measurements when a third 
neutralino can be measured, see section HHP I Now for less limited AA t ~ 200 GeV (which corresponds approx- 
imately to the range of variation we obtain once using all input, as we shall see later) those bounds are somewhat 
increased: 

300 GeV <h(Qewsb) < 410 GeV (6.6) 



Note that these bounds are very insensitive to the too variation within its accuracy range determined by squarks aS3 
sleptons in section 4. The limits ( 16.5b . (16.6b from EWSB constraints are anyway improved as compared to those 
obtained solely from two neutralino mass input in section 3. 



A. Naive tree-level counting of parameters 

Let us first start with a simple tree-level analysis in order to delineate the parameters entering in game. Concerning 
the lightest scalar Higgs mass, at tree-level, TO 2 '* ree is given by: 



ree = ]- \m\ + m| - ((mi + to|) 2 - 4mim| cos 2 2/3) 1/2 



Inverting this relation gives a (unique) solution for m A : 



m 2 h (m 2 z -m 2 h ) 
m? z cos 2 2(3 — to 2 



(6.7) 



(6.8) 



from which one can also derive: 



™A - (M 2 + m|/2)(tan 2 /? - 1) 
tan 2 (3+1 



m k = — — Vf ~ ^TiW " - — - > TO « d = 771 a - m k - V • (6-9) 



So very naively one may have thought to equate Eqs. ( 16.3b with Eq. ( 16.8b , deriving from this precise tan (3 constraints 
typically. But, as already mentioned, this tree-level analysis would be very unrealistic since the lightest Higgs mass rrih 
and m A get large radiative corrections, already at one-loop level, mainly from the stop/top loops, which are enhanced 
by a m\ dependence. Rather, as we did for the neutralino and squark sectors, we will try to incorporate the bulk of 
these large corrections consistently. 



B. Reconstruction of Higgs sector parameters in constrained MSSM 

The MSSM neutral Higgs mass-squared matrix reads, including radiative corrections: 

/ m| cos 2 (3 + vn\ sin 2 (3 + Su -(m| + m 2 A ) sin (3 cos (3+ S12 \ 
\ — (m| + m 2 A ) sin /3 cos /3 + S12 m| sin 2 (3 + m 2 A cos 2 (3 + S22 I 



(6.10) 



where Sij designate generically loop self-energy contributions. To explain our procedure, let us first consider a very 
simple approximation for nih (see e.g. ifsoll ): 



2, tree 



where 



3ffj m| 
S-K 2 m 2 w 



In 



Ml 



12M| 



(6.11) 



X t = At — [i cot (3 



M 2 



r 2 2 1 2 1 2 1 
= l m Q m t R + m t i m Q + 

i L , = [Ml - 4m 2 X 2 ]i 



(6.12) 



, ) + m 



(and 52 is the SU (2) gauge coupling). There is a large literature on various approximations to the MSSM Higgs 
mass(es) radiative corrections. We emphasize that in our actual analysis, we will not use the simple Eq. d6.1 U but 
a more elaborated expression, yet giving compact expressions for the needed Sij in Eq. ( 16. 1 Ob 1 391 . (This compact 
approximation is included in SuSpect as one possible Higgs mass calculation option). The latter incorporates some 



of the leading two-loop effects, but depends only on the very same parameters as in Eq. d6.1 U (apart from SM-Hke 
parameters). 

At this stage, if a rather precise information on the stop sector would be available, we could use an improved version 
of Eq. ( 16.81 ), incorporating radiative corrections, to derive ttia independently from the relation ( 16.31 ) 15 . However, 
within the SPSla input assumptions, this would not give any useful constraints (given also the poorly known tan j3 
constraints). We thus rather assume universality relations d6.4| i, such that Eq. ( 16.31 ), together with the m/, measurement, 
can determine both X t and tan (3. Since we also now assume complete universality of the squark sector, we can use the 
previous toq determination to calculate also the stop sector soft terms mi R and m^o (t he latter already obtained from 



the sbottom sector), as needed in Eq. ( 16.1 U (or its more complete generalization [39] actually used in the numerics 
below). The simplified expression of to/j we used is remarkably close to the full one-loop (plus leading two-loop) 
,001 ( 
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(in the DR scheme), with at most a 2 GeV discrepancy (and often much less) for a large range of 



MSSM parameters. (For the relevant SPSla case, one finds m 



simp 



111.28 GeV, while m 2 ' loop = 110.90 GeV.) 
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FIG. 19: Constraints on tan f3 and Xt = A t — /u/tan/3: indigo: one-a (68% C.L.) from a gaussian scan; green: same contour but 
for theoretical uncertainties Amj = 2 GeV instead. 



We now apply the previous analysis to the reconstruction of some of the parameters related to the Higgs sector for 
the SPSla test case with corresponding input mass accuracies. We perform again scan over parameters, taking into 
account the previous constraints on tuq available from squark and slepton sector (and letting X t arbitrary), deriving 
joint constraints on X t , tan (3. The variation of other parameters which it implies (like e.g. the variation of fi from 
EWSB relations ( 16.11 )) is consistently taken into account. This is illustrated first in Fig. |T9l where the range of points 
obtained from the prospected experimental accuracy: Am^(exp) ~ 0.25 from Table HIT1 are shown. Also shown on 
the figure is the similar range now corresponding rather to the present theoretical uncertainties in the Higgs mass 
calculations (we take the latter as Am/j(th) ~ 2 GeV IHOl 14711 ). Note again that the density of points in the tan /3, X t 
plane has no precise statistical meaning: the larger "density" seen for rather large tan (3 is an artifact of the logarithmic 
scale. Also the one-er band which is thiner for low tan/3 and small X t , is easily explained from the m/j dependence 



15 Note that even the full one- (two)-loop Higgs radiative corrections, i.e. with full expressions of the Sij above, do preserve a linear m? A solution 
in terms of similarly to Eq. J6.81 but with appropriate corrections from Sij. 
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FIG. 20: Constraints on tan j3 and Xt = At — [i tan (3: indigo: one-cr (68% C.L.); orange: additional points at two-cr (95% C.L.); 
green squares: remaining points (one-cr level) when adding sbottom mass measurements. 



on tan (3 and X t : more precisely for fixed tan 0, m,h(X t ) reaches a maximum for a certain X t (see e.g. Il40l0 . and it is 
more difficult to match the approximate SPSla value of mh for both small tan (3 and large (negative) Aq. The domain 
on tan (3 values obtained from this experimental Higgs mass measurement, at the one-cr level, is thus approximately 
determined as: 



5 < tan/3 < 40 . (6.13) 

These bounds are strongly correlated with Xt values, as the figure shows. For Xt one finds approximately: 

- 730 GeV <X t < - 450 GeV (6. 14) 

Note also that all the points shown with relatively large tan j3 ~ 35 — 40 are satisfying the previously mentioned 
theoretical constraint, with no tachyon fi , since when departing from the SPS 1 a benchmark, larger values of tan j3 can 
be compensated by values of Ao > —100 GeV. (This is of course not taking into account any additional experimental 
constraints on nth, m?, etc which, as already emphasized, are omitted here since they are more model-dependent 
constraints.) The resulting limits on X t and tan/3 are of course sensitive to rrih theoretical uncertainties as shown on 
the figures (though very little as concerns the tan j3 upper bound). In this case we obtain: 

4 < tan/3 < 40; -850 GeV < X t < - 330 GeV . (6.15) 

However the contours shown in Fig.[l9]did not take into account yet the additional constraints from the sbottom mass 
measurements which were obtained in previous section 5. This is now added in Fig. [20] (in green), where the one 
and two-cr domains obtained from mh measurement are compared. One can see that at the two-cr level the X t range 
increases slightly, for large tan (3, while tan (3 bounds are not much affected. On the other hand, the sbottom mass 
constraints put further upper limits on tan (3 < 27 — 28, in consistency with the previous analysis in section 5. 

These results are also compared with a top-down MINUIT minimization in Table IVIII using as input data rrih in 
addition to the gluino cascade sparticle masses (but not yet using the sbottom masses at this stage). One observes 
that the MIGRAD symmetric error on tan (3 appears to be very good, but comparing with the contours in Figs. [T9ll20l 
one can suspect that it is essentially influenced by the lowest limit on tan (3, and that the actual errors are particularly 
unsymmetrical due to the flat behaviour for large tan (3. This seems confirmed by the fact that the unsymmetrical 
MINOS positive error is not calculated by MINUIT. The lower tan (3 bound, on the other hand, is very consistent with 
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TABLE VII: Constraints on mSUGRA parameters mo, mi/ 2 ,Ao, tan j3 obtained from a standard top-down \ 2 fit of (first two gen- 
eration) squark, sleptons, neutralino/gluino masses, plus the lightest Higgs mass, for different choices of MINUIT minimizations. 



Data & fitted parameter 
(+ model assumptions) 


MIGRAD 
minimization (68%CL.) 


MINOS 
minimization (68%C.L.) 


nominal 
SPS la value 


m e - 2 ,m il ,+ 

m N 1 ^ m N 2 i m a +m h 

(1-loopRGE + noqR.C) 

q.l 

mi/2 
A 
tan f3(mz) 


(convergent) 
99.96 ±11.8 

250 ±4.9 
-100.5 ±150.0 

9.97 ±4.3 


(problems) 
99.9 t\o 5 
250 ±4.9 

q 07 + no limits 

-4.3 


100 
250 
-100 
10 



our alternative finding in Eq. ( 16.13l l. This behaviour is well explained by the flattest dependence for large tan j3, as 
clearly illustrated by the plots in Figs. [19] [20] For completeness we also show in Fig.[2T]the corresponding one- and 
two-(7 domains but in the (tan j3, A n ) plane, as well as the domain obtained once adding sbottom mass measurements. 
One can see that the sbottom mass bound have some impact on Aq bounds, at least for large — Aq, These results may 
be compared with the Aq limits obtained from the MINUIT fit in Table lVlfl 




1 10 20 , n 50 

tanp 

FIG. 21: Contours from a Gaussian scan in the plane (tan/3, Aq): indigo: one-er (68% C.L.); orange: two-cr (95% C.L.); green: 
remaining points (at one-a) when adding sbottom mass measurements. 



Thus the rrih measurement alone gives a mild upper bound on tan (3 in the case where only some knowledge on 
mAiQEWSB) is assumed (through sfermion-Higgs sector universality), but not on X t . The improvement on tan/3 
limits is mainly coming from adding the sbottom mass measurements. However this mild sensitivity, as well as the 
strong tan j3 — Xt correlation, is essentially determined by the variation of Aq which is assumed unknown at this 
stage. The uncertainties on fi (that are reasonable in this constrained MSSM case) have little influence (except perhaps 
for very small tan (3), since it is suppressed as /z/ tan/3 within X t . We shall see in the summarizing section |VID I that 
overall, the bounds are a bit tighter once combining all available informations, except those for Aq which depend quite 
much on the level of approximation used in RGE and radiative corrections. If an additional independent information 



on Aq would be available (e.g. if the stop masses could be measured) the constraints in the tan (3, X t plane could oi 
evidently improved. Better constraints on tan /3 are also prospected clearly if some of the heavy Higgses A, H, 
could be measured, that we do not assume in our 'minimal' set of input mass measurements and within the SPSla 
benchmark case. (It would be rather straightforward, however, to extend our above analysis in case such additional 
information would be available.) 

We finally mention that our results are not drastically changing once using the more complete two-loop level ra- 
diative corrections for mh (assuming in such case universality/mSUGRA relations to calculate all higher order cor- 
rections), except for the bounds on Aq that are substantially worsened in this case, as will be confirmed by further 
MINUIT fit results summarized in section [VI Dl This comes from an observed flattest dependence of nih(Ao) for 
large — Aq at the two-loop level, in the range considered for the other relevant parameters. As far as we can see, this 
may be explained by the fact that within the two-loop level radiative corrections to m/J^, 51], there is somewhat 
much room for possible cancellations of the Aq dependence entering in higher orders, most notably in all squark and 
stau contributions both at the one- and two-loop level. This is another example of the caution needed to interpret 
results when performing minimization constraints using different level of theoretical approximations. 



D. Combining all information and comparison with standard fits 



TABLE VIII: Combined best constraints from all gluino cascade decay sparticles and lightest Higgs mass measurements for the 
bottom-up approach in different MSSM scenarios. (NB a star in superscript indicates discrete ambiguities in the reconstruction in 
the non-universal MSSM case, as discussed e.g. in sub-sec. 3.1 and Appendix A). 



Model or assumptions 


Parameter 


Constraint (GeV) 


nominal SPSla value 
(from SuSpect 2.41) 


general MSSM 


Mi(Qewsb)* 


~95-115 


101.5 




M 2 {Qewsb)* 


~ 175-225 


191.6 




M 3 (Qewsb) 


~5 80-595 


586.6 




(l™l L + ^) 1/2 (Qgut) 


~68-89 


(|) 1/2 100 ~ 79 




m Q3L {Qewsb)* 


~488-518 


497 




m,b R (QEwsB)* 


~5 10-540 


522 




h{Qewsb)* 


~280-750 


357 




t&n{3(Q E wsB) 


~l-36 (th. bounds) 


9.74 


+ m N i 


h{Qewsb)* 


~350-372 


357 




tan P(Qewsb) 


~2.7-36 (th. bounds) 


9.74 


q, /-universality 


ml' l {QGUT) 


~90-112 


100 




m Q3 L (Qewsb) 


-490-506 


497 




m bR (Q E wsB) 


~5 13-530 


522 


Mi,i = 1, ..3-universality 


Mi(Q G uT) 


~ 245-255 


250 


bi , &2 +universality 


tan P(Qewsb) 


~3-28 


100 


mSUGRA 


mo 


—90-112 


100 




mx/2 


~245-255 


250 




-x t 


~450-730 


530 




-Aq 


~ -100-350 


100 




tan f)(mz) 


~ 5.5-28 


10 



We finally summarize and combine all previous constraints from the different sectors, both in general MSSM or 
with universality relations, in Table I Villi By "combined best constraints" we simply mean to evaluate the joint limits 



by crudely superposing the constraints obtained in the three different sectors of gaugino/Higgsino, squarks/slepto'ril 
and Higgs, which eventually results in slight improvements on some parameter limits. We did not attempt to perform 
a more elaborated statistical combination of the three sector constraints within our present approach. However, these 
results are compared with those obtained from a top-down x 2 fit of basic mSUGRA parameters to the same mass 
measurements in Table HXl As one can see, the results are overall qualitatively consistent, except perhaps for the Aq 
limits. The fact that Aq bounds (and to some extent tan/3 bounds as well) are worsen when fitting parameters in 
mSUGRA at the full two-loop level instead of one-loop, was explained previously due to the possible cancellations 
of the Aq dependence in higher order contributions. Note that for tan /3, the MIGRAD symmetric minimization 
error found here appears very optimistic as compared to the upper bound derived from the bottom-up reconstruction, 
tan/3 < 28. As already mentioned above it is possibly much influenced by the lower bound that results from the 
combination with the Higgs sector measurements. The MINOS upper bound is however more consistent with the 
bottom-up result. 

TABLE IX: Combined constraints on mSUGRA basic parameters from a standard top-down x 2 fit with MINUIT of all gluino 
cascade decay sparticle masses plus mj, measurements. 



Model and assumptions 


Parameter 


68% C.L. limits (GeV) 


SPS la value 


mSUGRA 


mo 


99.96 ± 11.2 


100 


2-loop RGE + full q R.C. + 2-loop m h 




(99.95 ±11.7) 




(1-loop RGE+ no q R.C. +simple m% R.C.) 


m 1 / 2 


250.0 ± 3.7 
(249.5 ± 4.7) 


250 




A 


-104.2 ±379 
(-100.6 ±136) 


-100 




tan/3(mz) 


9.9 tl A 7 
(9.96 ±4.11) 


10 



VII. CONCLUSION 
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We have examined some specific bottom-up reconstruction strategies at the LHC, both for general and universality- 
constrained MSSM parameters, starting from a plausible set of incomplete measurements of a few MSSM sparticles. 
Using sparticle mass measurements mainly from cascade decays of gluino and squarks, and the lightest Higgs boson 
mass, we have constructed different algorithms, based on rather simple semi-analytical inverse relations between the 
MSSM basic parameters and mass spectrum, incorporating radiative corrections, when known, at a realistic level. 
We have determined constraints on the relevant basic MSSM parameters from the expected mass accuracies, under 
different theoretical assumptions on the degree of universality of some of the parameters. We have also exhibited 
analytically the possible discrete ambiguities in the reconstruction of some of the basic parameters, when using only 
mass input in a general unconstrained MSSM without extra knowledge or assumptions on the relative hierarchy of the 
relevant parameters. This is the case in particular for the gaugino/Higgsino sector parameters. These constraints have 
been also compared at different stages in a sector-by-sector analysis, with those obtained from more conventional 
top-down approaches of fits to data with minimization procedures. The results are overall consistent, which is an a 
posteriori check that our rather naive semi-analytic approach, with many approximations, does essentially capture the 
sensitivity on parameters. 



Regarding the SPSla reconstruction example studied here, more quantitatively we have shown that a rather 
limited data set, consisting of merely the measurement of sparticle masses involved in gluino cascade decay to a 
few percent accuracy, may still provide reasonably good constraints on some of the relevant MSSM parameters. 
This is in particular the case for the gaugino mass parameters and the squark and slepton soft mass terms, even 
for an unconstrained MSSM without universality assumptions. If a precise measurement of the lightest Higgs 
mass is available, additional constraints (though moderate) on tan/3, X t = A t — [ij tan/3 are obtained, but only if 
Higgs-sfermion mass term universality is assumed. Interesting constraints for a non-universal MSSM are however 
more challenging to obtain in general for the Higgs sector parameters (as well as for tan (3), unless more precise 
measurements would be available (or data from another sector, like heavy Higgses and/or scalar top masses typically), 
as is also known from other analysis]^. We stress again that considering only the gluino/squark cascade (plus the 
lightest Higgs) sparticle identification in our analysis is not motivated by a strong prejudice against other potential 
SUSY-discovery processes at the LHC. Indeed the cascade in Eq. ( 12.11 ) may be considered already quite specific from 
a general MSSM viewpoint, but it gives us a well-defined and rather minimal input set for testing our approach, 
and comparing it with other analyses for the very much studied SPSla benchmark. Even if our results are probably 
not very new to the experts, they illustrate that a step-by-step semi-analytic approach can help to exhibit better the 
sensitivity of basic MSSM parameters with respect to given sparticle mass or other data, which may be more difficult 
to grasp from global top-down fits. In many cases, the non-linear and non-symmetric behaviour of error propagations 
is exhibited by our analytic bottom-up approach in a more explicit way (as illustrated typically for tan (3, comparing 
Tables IVIIII and IIXI) . and similarly for the possible discrete ambiguities e.g. in t he g augino/Higgsino sector. As 



compared to other recent analysis of MSSM constraints at LHC (e.g. |5|, [a [17|, 12511 ). our results are difficult to 
compare quantitatively in very detail, since the data used are often different (with generally more input sparticle 
masses assumed in most other analysis). We find however a rough consistency on the expected sensitivity of the basic 
MSSM parameters, as discussed above. 



The relatively simple algorithms described here are rather flexible, and may be easily interfaced with more elab- 
orated simulation tools, that would allow in particular a more realistic statistical treatment of the different sources 
of experimental and theoretical uncertainties. Moreover, some of our analytical relations may be used at least as 
"Bayesian priors" guideline to other analysis, in a way complementary to the Markov chain techniques. We plan 



indeed a more refined statistical analysis in the future |)41|], by possibly combining our approach with Bayesian ana 
Markov chain techniques of ref. [TtL Eill . 

Our approach may thus provide a useful complementarity to more elaborated simulations, as well as possible cross- 
checks. Moreover it is not at all restricted to the LHC phenomenology: some of the algorithms described here, for 
instance in the neutralino and squark/slepton sectors, may be readily used for ILC data upon straightforward changes 
in sparticle mass accuracies. However, a similar approach for the ILC, following e.g. the studies made in refs. 1 19112311 . 
deserves specific analysis beyond the scope of the present paper, due to the different sparticle spectrum expected to be 
reached at the ILC, which will imply slightly different inversion algorithms. 
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APPENDIX A: SOLVING THE GAUGINO/HIGGSINO PARAMETERS FOR DIFFERENT INPUT 
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In this appendix, we give for completeness explicit analytic solutions relating neutralino masses and gaug- 
ino/higgsino parameters, depending on different input and output choice and depending on general MSSM or more 
constrained universality assumptions, relevant to sections 3.1-3.4. 

1. Two neutralino mass input 

From Eqs. ( 13.3b , ( 13.4b , one can solve Mi , M2 for input masses , m^, plus fi, tan (3 input: after expressing Mi 
e.g. from Eq. ( 13.3b (see also lfl8ll ): 



P12 [M 2 + m% + P12 + (M 2 - S 12 )S 12 ] + M 2 11 m 2 z s 2 w sin 2/3 
S 12 P 12 - M 2 (^ + P 12 ) + c? w m%»sm2P 



(with S12 = m ] y i + mjY 2 , Pi 2 = in^^mfj^), M 2 is given by a standard second order equation with solution: 



M 2 = ^ ± ^ - 4aiCl (A2) 
2ai 

where +^777 corresponds to Mi < M 2 and —^/~~- to M 2 < Mi respectively, and 

ai = P12 [m 4 + Pi 2 (Pi2 + m%Sw) + M 2 (2Pi2 - Sf 2 + m 2 z s 2 w ) + ^m 2 z s 2 w sm2[3 S 12 ] 
h = Pi 2 [S 12 {-^ + ^ 2 {~m 2 z -2P 12 + Sl 2 )-P 12 {P 12 +m 2 z {s 2 w -c 2 w )) 

-/xm|(Pi2 + {S 2 12 - Pn){sw -4v) + 2n 2 c 2 w )sm2f3] 
ci = Pi 2 (^i 2 (M 2 + m% + Pi 2 )(^ 2 + c 2 w m l + A2) - (^ 2 + Cwm 2 z )P 12 Sl 2 

+lim 2 z sin2/3(S*i 2 (Pi2 + ^(m 2 + m| + 2P i2 - S 2 l2 )) + c 2 w m 2 z fx sin 2/3)) . (A3) 

The occurrence of two M2 solutions in Eq. iA2i is due to the structure of Eqs. 0.31 >. ( 13.4b . but this ambiguity cannot 
be resolved from the sole knowledge of two neutralino mass input, since all the available relations from Eqs. (13.2b are 
completely symmetrical under m# permutations. Typically in the simple limit of neglecting all Z?-terms i.e. 

mz — » in the neutralino mass matrix Eq. ( 13.1b , the four neutralino mass eigenvalues are trivially given as 

m# 4 (i = l,..,4) = Mi,M 2 ,| / u|,| j u| . (A4) 

with all possible permutations i.e no ordering implied. The exact mass eigenvalues with mz 7^ are no longer 
invariant under e.g. Mi «-> M 2 permutations, but concerning the inverted relations, even if all physical masses were 
known and thus ordered according to e.g rrif Ji < ■ ■ ■ < m^, there is no way to determine from the sole knowledge 
of a given mass if it corresponds to a Bino, Wino or Higgsino. This needs information on diagonalization matrix 
elements, that can be partly accessed if measuring the couplings of neutralinos to other particles. (However, for known 
/i input, measuring a third neutralino mass could in principle solve the Mi , M2 ambiguity, since the two different 
solutions gives different m^ 3 , m^ 4 (which are most simply obtained by solving Eq. ( 13.61 ) as a second order equation 
for e.g. m^J. But in practice this needs an accurate knowledge of both [i and the three neutralino masses). 

In Fig.|22]the two Mi, M2 solutions are illustrated as functions of fi (see respectively M 2 (l, 2) and Mi(l, 2)), with 
the two neutralino mass input and tan/3 fixed to SPSla values for reference. Note that in Eq. ( lA2b there can be values 
of /1, tan/3 for which ai — » 0: for example for SPSla values of tan/3, ai(/x, tan/3) has four zeros for real values 
of /1: [i ~ —165, —110, 118, 157 GeV. However these poles with respect to [i are artifacts since for a\ = strictly, 
the (unique) solution reduces to M 2 = —c\jb\ (moreover the possible zeros of bi(fi) do not coincide with those of 



ai (//)). Nevertheless sufficiently close to these [i "pole" values, one of the two Mi solutions can become arbitrarily 
large, as can be seen on Fig. [22] 16 . More physically, this means that such large values of M2, with their corresponding 
/i pole values, are still consistent with the same neutralino masses , . 

As such, Eq. dA2b is not very illuminating: it may be useful to expand to first order in the D-term, i.e. to 0(m 2 z ), 
which gives a very simple expression: 

„ „ (mcr + u sin 2/3) 
Ah = m# a + c 2 w m 2 z N > - 2 - (A5) 

forthecaseMi < M% (and the same expression with m# — > for the case M2 < Mi). (For M\ one better keeps 
the exact and relatively simple expression Eq. (IA1I) ). Here the pole at fi = in Eq. dA5t is clearly an artifact of 
first order in m 2 z expansion, and simply reflects that this approximation is not valid for /x ~ m^.. Eq. d A51 > is certainly 
a good approximation for /j, ^> , m^ 2 and m^, ^> mz- Moreover it can still be a rather good approximation 
even for moderate \i and neutralino mass values. For instance for SPSla values in Table HH and neglecting self-energy 
radiative corrections to neutralino masses one has fj, ~ 357 GeV, ~ 176.2 GeV and M| a:act (SPSla) ~ 191.6 
while M2 Pprox (SPSla) ~ 192.6 which is only .5% difference, i.e. less than the magnitude of radiative corrections to 
neutralino masses, which are about ~ 4 GeV for m ^ 2 . We have used however the exact expressions iA2i throughout 
our analysis in section 3. 




FIG. 22: The two Mi, M2 solutions (respectively M^(l, 2) and Mi(l, 2)) from Eqs. dA lb . iA2i for two neutralino mass rn^ i , 
rriff 2 input, as functions of /1. The dots and triangles corresponds to the solutions once a third neutralino mass m^) input is 
used, completely determining (for fixed tan f3). The neutralino mass input values and tan j3 are fixed to their SPS la values for 
illustration. 



In fact only two of the ai (fi) zeroes lead to this pole behaviour of M2(fJ.), as seen on Fig. l22l (e.g. at fj, ~ —165 and 157 GeV for M2 > Mi 
and the two other zeros for M2 < Mi), whereas the other zeros of ai give bi((j.) > such that the correct behaviour M2 — » — C1/61 for 
01 — > is recovered. 



2. Three neutralino mass input: general case 
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Knowing three neutralino masses allow to use one more equation from d3.2b to solve e.g. for M±, M2 and /i in terms 
of the three neutralino masses and tan (3 input values. (Alternatively, we could also solve this system to determine 
rather tan (3 for fixed fi, but since it only depends on sin 2/3, any tan (3 dependence drops out for large tan {3. Thus we 
can anticipate without calculations that there can be no interesting tan (3 upper bounds, irrespectively of the number 
of neutralino masses measured, given their expected LHC accuracies). The extra relation to determine /i 2 is Eq. ( 13.6b 
obtained from the trace and second invariant in Eq. ( 13.21 . Note however that the sign of fj, is still not determined. 
Combining Eq. ( 13.61 ) with Eqs. d3.31 l. ( 13.41 ) gives a high (sixth) order polynomial equation in /i (or equivalently in M\ 
or M2) which cannot be solved fully analytically. It is however easy to solve iteratively using e.g. Eq. ( 13.6b on the 
solutions dA2b . dAll ) (upon having chosen a definite order for the relative values of M\ and M2). This iterative solution 
converges very quickly. The solutions for Mi, M2 and [i obtained for tan/3 input fixed to SPSla value are illustrated 
by the different dots and triangles in Fig. [22] 

Alternatively it is useful to consider, as in the case of two neutralino masses above, approximate solutions obtained 
by expanding to first order in m 2 z . This gives very simple expressions: 



Mi = m 
M 2 = 



( m Ni + 111 n 3 sin 2/3) 



Nl T- o w ,„ z 3 _ 2 

N 3 Ni 



(m N2 + m Na sin 2/3) 



-w" l z 2 2 

N 3 N 2 



\lA = m^--^(l + sin2/3 L_p ^ ™ N * (A6) 

2 ( m N 3 ^ m N 1 > ( m N 3 - m N 2 ) 

corresponding to the first case Mi < M2 (i.e. a Bino LSP), and similar expressions with <-» m N for the case 
M2 < Mi. Actually, considering three neutralino mass input and Mi,M 2 ,/z output adds more discrete ambigui- 
ties than in the two neutralino mass case: accordingly in a most general MSSM, without further knowledge on e.g. 
neutralino couplings, or further theoretical assumptions, one should in principle consider all possible ordering among 
Mi, M2, fj, values, i.e. six possible cases (where the LSP mass m Ni can be either Bino, Wino or Higgsino). It is 
straightforward to derive such other solutions by appropriate permutations of the three m N ., within Eqs. dA6K For 
all cases Eqs. (IA6b are very good approximations, at least as far as are not small compared to mz, e.g. for the 
SPSla case typically the differences with exact results is of order or below the percent level. 

3. gaugino mass universality: two or three neutralino mass input 

Assuming gaugino mass universality one can solve again Eqs. ( 13.3b . d3.4l ) but changing input/output: now the EWSB 
scale values of Mi and A/2 are given from M3 by Eq. ( 13.5b . so that one can determine and tan (3 assuming two 
neutralino mass input from a linear system, which gives 

/i 2 = -— (A7) 

a-2 



with 



P (M J- JUT c 1 A2-M1M2 , 

a 2 = P\2 (Mi + M 2 - S12 H — ) 

M\2 

b 2 = P12 [P12 (Mi + M 2 -6-12) + 5 12 ((Mi - m^)(M 2 - m^) - (M x + M 2 - m^)m 



2\, J L ^12 i J 12 + (AJ 1 -Sl 2 )(Sl2-M 2 ) . 
~ JV2 - H + ^ ( ^ + I 




A ! 2 
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. [P? - M 1 M 2f x 2 + P12 ((Mi + M 2 )S 12 - MiM 2 - Sf 2 + fi 2 + m|)] 

sin 2/3 = — t 2— i- (A9) 

/x mr z M\ 2 

with Mi 2 = CyyMi + s 2 v M 2 and other notations as in Eqs. (|Alb - (IA3t . Note that we assume \i to be real. The case of 
three neutralino mass input is treated by combining Eq. (I3.6b -Eq. 



APPENDIX B: BOTTOM-UP RENORMALIZATION GROUP EVOLUTION IN SUSPECT CODE 

In this section we illustrate in some details an important ingredient of our bottom-up reconstruction procedure 



which is the renormalization group evolution (RGE) 12a 12711 of MSSM parameters from low to high energy. We take 



this opportunity to present the results from an available option of the SuSpect code, illustrating some properties of 
this bottom-up RG evolution which are quite general and independent of the present analysis and LHC cascade de- 
cay phenomenology. (This option has recently been adapted to the more suitable "Les Houches accord" input file 



conventions |52J.) 

The fact that the RGE for the MSSM parameters are "invertible", i.e can be evolved from high to low scale and 
backward, is a rather obvious feature of any such coupled differential equations. However, what makes it less straight- 
forward in the standard approach to MSSM calculation with RGE, is that there are actually three (at least) different 
energy scales in the game, with corresponding boundary conditions: 

• the high (generally GUT) scale, at which one defines e.g. the mSUGRA parameters with eventually some 
universality relations; 

• the (low) electroweak symmetry breaking scale, at which the soft breaking and other relevant MSSM parameters 
are evolved to in a standard top-down evolution; 

• finally the scale n%z, or eventually some other low energy scale, which enter as other boundary conditions e.g; 
via the precise measurements of the gauge couplings. 

The i nterplay between these different scales in MSSM spectrum calculation codes (such as SuSpect J3II1 and similar 



codes! 28. 129L 13010 needs, among other things, iterative algorithms for the evolution between the three different scale, 



with a consistent implementation of any possibly known radiative corrections to the parameters (gauge couplings, top, 



bottom masses, etc) entering as boundary conditions to the RG evolution. (We refer to the SuSpect manual 13111 for 
details on such RGE algorithms.) 

Such a bottom-up evolution option has been available in a beta-version for some time in SuSpect but was not 
publicly available nor illustrated until now. In Table IXl we illustrate for the SPSla benchmark point this bottom-up 
RG evolution of all parameters. We show in particular some important features on the error propagation in such a 
procedure, and which parameters are more sensitive to this dispersion. The input parameters (in the second column) 
were obtained in a first stage from a standard (top-down) run from the mSUGRA SPSla input parameters in Eq. (12.2b 
(with m top = 175 GeV). The third column gives the corresponding output values resulting from a RG evolution up to 
a GUT scale from a bottom-up SuSpect run under the most general MSSM assumptions, i.e. without any a priori on 
possible universality relations at the GUT scale. One can already notice from this that the agreement with the original 
mSUGRA parameter values is excellent: the discrepancies are of order 0(1O~ 3 ) that are in fact consistent with the 
accuracy chosen (i.e. the intrinsic error of the numerical evolution of the RGE as performed with a Runge-Kutta 



algorithm, as well as the intrinsic error due to necessary iterative procedures 113 111 ). We point out however that the 
results shown here correspond to the choice of m,H u , m,H d input, while the correspond results for rriA, M input choice 
are a little bit worse for to# u (Qgut) an d m H d (Qgut), which is attributed to a certain precision loss in our algorithm 



when passing from one to the other input set, that involves calculations and iterations inclu ding two-loop radiative 1 
corrections. (This is in particular due to the well-known Tr[YM 2 ] terms present in the RGE 026L 12711 . which involve 
a combination of soft scalar terms, including m? H and mf^ : this combination vanishes by definition in mSUGRA, 
and remains zero at all scales in a top-down RG evolution. However, when using bottom-up RGE the approximate 
m 2 Hu , rrijj values as obtained from tjia, £t from EWSB conditions Eqs. ( 16. Il l does not satisfy Tr[YAI 2 ] = exactly, 
and this induces a not completely negligible departure when evolving from low to GUT scales.) 

Next, the fourth and subsequent columns give the deviations in the output parameter values corresponding to ±1% 
deviations of some relevant parameters input values: M3, mjj„, ttiq3 l , that we chose on purpose as they give the 
most important sensitivity to other parameters deviations. One can see that the deviations induced in other parameters 
remain generally reasonable, at the percent level, notably for the gaugino masses, and most of the squark and sleptons 
(except for the one sfermion mass which is varied in each case). In contrast, the deviations induced on the scalar mass 
m,H u can be huge, even for a moderate percent deviation, i.e. there is a large amplification or dispersion of error. 
This is of course explained from the detailed RGE dependence of rriH u on other soft terms, resulting in a very strong 
sensitivity on certain other parameters. The same is true to some extent for m tR . The strong sensitivity of mjy u e.g. 
on the top mass value and other parameters through its RGE is a well-know feature, but more precisely our illustration 
here indicate that it will be very difficult to have a precise determination of to# u at a high GUT scale, even in the 
very optimistic case where all other MSSM parameters would be know quite precisely. Consequently it will be very 
difficult to check for eventual universality of the soft scalar mass terms with the squark and slepton soft terms. One can 
however turn this argument the other way round, and deduce from this illustration that, while to# u plays a crucial role 
in the (radiative) electroweak symmetry breaking, the sparticle masses at low scale determining most of the collider 
phenomenology are not very much dependent on its precise value at GUT scale. (This is somewhat similar here to the 
"focus point" properties observed for large too values in other mSUGRA parameter choices). More generally, many 
features as those illustrated in Table IXl may be important to keep in mind for any realistic bottom-up procedure. The 
above bottom-up RGE procedure had been used in several stages of our analysis, as indicated in the main text. 
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TABLE X: Bottom-up RG evolution of SPSla parameters from SuSpect 2.41 with illustration of error propagations. Input pa- 
rameters (2nd column) were obtained from a standard (top-down) SuSpect run from SPSla input parameters in Eq. l !2.2t (with 
m t op = 175 GeV). The third column gives corresponding output values once evolved back to a GUT scale from a bottom-up 
SuSpect run with general MSSM option. Fourth and subsequent columns give deviations in the output parameter values corre- 
sponding to ±1% deviations of some relevant parameters input values as indicated. (Variation range is given explicitly when 
non-symmetrical) 



par. 


input(GeV) 


GUT output 


AM 3 = q=l% 


Am Hu = =Fl% 


AmQ3 t = =Fl% 


Qewsb 


465.5 


~ 2.47 10 16 


0.1 % 


0.1% 


0.1% 


Mi 


101.5 


250.004 


negl. 


negl. 


negl. 


M 2 


191.6 


249.998 








M 3 


586.6 


249.999 


±2.2 






m H d 


(179.9) 2 


(100.004) 2 


(100.6) 2 - 


(100.7) 2 - 


(101. 2) 2 - 








(99.4) 2 


(99.2) 2 


(98.7) 2 




-(358.1) 2 


(100.017) 2 


(132.6) 2 - 


(64.9) 2 - 


(63.7) 2 - 








(48.4) 2 


(124.4) 2 


(126.4) 2 


(m P A° lB ) 


398.8 










(ft) 


356.9 


353 








rrie L 


195.5 


100.004 


100.2-99.8 


100.8-99.2 


101.5-98.5 


m TL 


194.7 


100.004 


100.2-99.8 


100.8-99.2 


101.5-98.5 


m eR 


136 


99.998 


100-99.9 


98.4-101.6 


96.8-103.1 


m TR 


133.5 


99.998 


100-99.9 


98.4-101.6 


96.8-103.1 


mQi L 


545.8 


100.001 


121-72 


99.7-100.3 


99.1-100.8 


m Q3 L 


497 


100.005 


131-52 


94.6-104.6 


55.2-130.4 


m UR 


527.8 


99.997 


121-72 


101-99 


101.8-98.1 


m tR 


421.5 


100.006 


140-14 


90.6-107.5 


81.9-115.3 


m dR 


525.7 


99.997 


121-72 


99.4-100.6 


98.7-101.3 


m bR 


522.4 


99.997 


122-72 


99.4-100.6 


98.5-101.5 


-A t 


494.5 


100.009 


111 - -89 






-A b 


795.2 


100.002 


106 - -94 






—A T 


251.7 


100.002 


100 - -99.9 








677.3 


100.005 


108 - -92 


negl. 


negl. 


-A d 


859.4 


100.001 


105 - -95 






-A e 


253.4 


100.002 


100 - -99.9 






tan (3 


9.74 
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